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Chapter 1 


Introduction 
1.1 Background 

Noise is an inevitable conseqnence of rotary wing flight, and can severely 
constrain helicopter operations. Civil helicopter operations are frequently limited by 
the public acceptance of aircraft noise, as well as local and national regulations on 
permissible noise levels. These noise criteria restrict the use of helicopters for pas- 
senger transport, sightseeing, electronic news gathering, and aerial law enforcement. 
Noise criteria also restrict the construction of heliports, preventing the adoption 
of point-to-point vertical lift commuter aircraft, often touted as a solution to the 
overcrowding of airspace worldwide. Military operations are likewise adversely af- 
fected by rotorcraft noise. While rotorcraft frequently evade detection by visual, 
infrared, or radar means by flying nap-of-the-earth, low to the ground and amongst 
terrain features, the lower frequency acoustic waves produced by rotorcraft carry 
over long distances and can refract around terrain. Consequently, rotorcraft are 
often hrst detected by acoustic means. Acoustic detection can be achieved using 
anything from elaborate arrays of electronic detectors to the unaugmented human 
ear, and can greatly decrease the survivability of helicopter missions. In both civil 
and military cases, the acoustic impact of helicopter operations can be reduced 
through either the design or the operation of the rotorcraft. However, changing the 
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design of the helicopter to reduce noise typically results in a reduction in other key 
performance parameters, such as payload or maximum speed. In addition, design 
changes offer no remedy for rotorcraft which have already been produced and are 
currently in widespread service. Designing helicopter operations for reduced acoustic 
impact is a realistic and promising avenue for increasing the number and variety of 
helicopter operations without exceeding existing limits on community annoyance 
or survivability. 00 Designing helicopter operations for reduced noise requires an 
accurate method of assessing the acoustic impact of individual helicopter operations. 


1.2 Ground Noise Contours 


The acoustic impact of helicopter operations is typically quantihed using ground 


noise contours. (Figure l.f on the following page ) For community noise, ground 


noise contours are used in order to plan helicopter operations so as to minimize 
annoyance in noise sensitive areas, to ensure compliance with community noise 
regulations, and to plan the location and design of heliports. Ground noise contours 
are also an increasingly valuable tool for military mission planning — designing flight 
paths that minimize the chance of aural or electronic detection, thereby improving 
the survivability of the helicopter. 

The process of generating ground noise contours requires three models; a 
helicopter noise source model, an atmospheric propagation model, and a receiver 


model. (Figure 1.2 on page 4 ) The helicopter noise source model describes the 


external noise radiation of the helicopter as it moves along the desired flight trajectory. 
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Figure 1.1: A notional ground noise contour for a helicopter on an 
approach trajectory. 


The atmospheric propagation model describes how the noise radiated by the helicopter 
is transmitted from the helicopter through the atmosphere and to observers on the 
ground. The receiver model describes how the noise heard on the ground impacts the 
observer. Atmospheric propagation models have received much attention recently 
with high-hdelity physics-based models being incorporated into the current generation 
of ground noise contour estimation tools. |1] While most current aural receiver models 
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are based on empirically generated descriptions of the freqnency response of the 
human auditory system, [Hj some recent research has suggested that these models 
may be inadequate for harmonic noise sources and that time-domain models, which 
capture both the magnitude and phase of the noise as a function of frequency, may 
be required. [B] The focus of this dissertation is on helicopter noise source modeling: 
accurate noise contour estimation is only possible with an accurate model of the 
noise at the source, no matter the quality and sophistication of the propagation and 
receiver models. 



Figure 1.2; The three components of ground noise contour modeling. 
From left to right; the helicopter source noise model, the atmospheric 
propagation model, and the receiver model. 


1.3 Helicopter Noise Source Modeling 

Helicopter noise source modeling can be broadly divided into two approaches, 
theoretical models that develop estimates of helicopter external noise radiation from 
hrst principles and empirical models which incorporate measurements of helicopter 
noise in the construction of the model. The state-of-the-art in theoretical modeling 
is still very time consuming and requires validation for new conhgurations and 
operating conditions; routine theoretical modeling of helicopter noise for specihc 
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configurations over a wide range of operating conditions is not yet practical. On the 
other hand, numerous empirical modeling schemes are currently in use. The two 
types of empirical models in general use are integrated models and source hemisphere 
models. 

1.3.1 Integrated Noise Models 

Integrated models have been used in the Federal Aviation Administrations 
Heliport Noise Model (HNM),[7] Integrated Noise Model (INM),[8] and the upcoming 
Aviation Environmental Design Tool (AEDT).|9] Integrated noise models are devel- 
oped from measurements of noise on the ground captured during aircraft flyovers. 
Integrated models provide a representation of the time-averaged noise radiated to 
points on the ground for an aircraft flying along a theoretically inhnite flight track. 
The variations of the integrated noise levels with respect to the distance between 
the source and receiver are described by empirically developed Noise-Power-Distance 
(NPD) curves. [TOj NPD curves are developed for three standardized flight conditions, 
also used for noise certification: level flight, takeoff, and approach. The distribution 
of noise levels by frequency is handled by assigning a spectral class to the NPD 
curves where the “spectral class” describes a normalized frequency spectrum that 
is assigned to the noise radiated by the noise source model in all directions for all 
flight conditions. 

In an extension to the modeling method for helicopter noise, some measure 
of the lateral variation in noise levels is achieved through the use of separate NPD 
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Figure 1.3: A diagram of the Integrated Noise Model three microphone 
measurement standard. 


curves for the left and right sides of the helicopter. These NPD curves are developed 
from measurements of the helicopter using three above-ground microphone locations 


for each of the three standardized flight conditions, as shown in Figure 1.3 However, 


integrated modeling does not explicitly model the overall directivity of the noise 
radiated from the helicopter. Lateral directivity is coarsely modeled by differentiation 
in the left, right, and center NPD curves as well as an empirically determined sideline 
correction, which extrapolates the left or right NPD values for observers far from the 
flight path of the helicopter. Longitudinal directivity is not directly captured by the 
model, but is implicit in the distance dependence of the NPD. Similarly, variation 
in noise levels with changes in flight condition is coarsely modeled using data from 
only the three standardized noise certihcation flight conditions. Integrated modeling 
is suitable for estimating averaged noise levels over many repeated flights — and is 
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considered applicable over only a limited geographical area, such as an airport — but 
does not accurately capture the noise radiation from a single flyover event. Integrated 
Noise Model was initially developed for fixed-wing aircraft, where due to the large 
volume of air traffic, time averaged noise levels provide a useful measure of the 
noise impact of airport operations. However, due to the typically lower volume of 
rotorcraft flight operations, and the unique character, directionality, and sensitivity 
to operating condition of rotor noise sources, integrated noise modeling may not be 
well suited to assessing the noise impact of most rotorcraft operations. 


1.3.2 Source Noise Hemisphere Modeling 

In order to address these deficiencies in modeling the external noise radiation 
of rotorcraft, several more sophisticated simulation methods have been developed 
based on the concept of empirically constructed source noise hemispheres (Figure [L4| 


on the next page). The fundamental assumption of these methods is that for far- 


held observers, the helicopter can be represented by a single compact noise source 
occupying a point in space, typically chosen to be the main rotor hub of the helicopter. 
As the helicopter hies over a ground-based array of microphones, the microphones 
measure the noise along the set of directions between the assumed source and 
the microphones. The noise can then be “de-propagated” along these directions, 
accounting for the time-delay between the emission of noise from the assumed source 
and the reception of noise at the microphone, to points a hxed distance away from 
the compact source by accounting for transmission losses due to spherical spreading 
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and atmospheric absorption. The result is a measurement of the noise levels on the 
surface of a hemisphere set a hxed distance away from the main rotor hub of the 
helicopter. Typically, the noise levels on the surface of the hemisphere are given in 
terms of 1 /3-octave or narrowband frequency spectra. This provides a description of 
the noise radiated by the helicopter for a single measured flight condition in terms 
of amplitude, frequency and direction. By repeating the process for various flight 
conditions, a database of acoustic hemispheres can be generated for use in ground 
noise contour generation. 



Figure 1.4: An source noise hemisphere centered about the main rotor hub. 


1.3. 2.1 Rotorcraft Noise Model 


The most popular hemisphere-based helicopter noise-modeling tool in the 
United States is the Rotorcraft Noise Model (RNM)[TT], developed by NASA, 
DoD, and Wyle Laboratories. RNM develops noise hemispheres from acoustic 
measurements of a helicopter in steady straight-line flight over a linear array of 


microphones, as shown in Figure 1.5 on the following page As the helicopter passes 
over the array, each microphone captures a range of directivity angles over the 


duration of the measurement (Figure 1.6 on page 11). Since the acoustic state 


of the helicopter is assumed steady, these measurements are used to construct a 
single hemisphere representing one flight condition. Flight conditions are defined in 
RNM in terms of the airspeed of the helicopter and the flight path angle. In order 
to characterize the helicopter, numerous combinations of airspeed and flight path 
angle are tested, resulting in a large database of acoustic hemispheres for steady 


straight-line flight conditions (Figure 1.7 on page 12). These hemispheres may then 


be used to estimate noise levels on the ground by dividing the desired trajectory into 
a sequence of steady straight-line segments. Each segment is then associated with a 
measured hemisphere — if an exact match for a desired flight segment is not found 
in the database of hemispheres, a new hemisphere is generated by interpolating the 
noise levels on the surface of the hemisphere between known values of airspeed and 
flight path angle. For turning flight segments, RNM has the capability to rotate the 
noise source hemispheres by the bank angle of the helicopter, but makes no other 
adjustments to the noise source. | 
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Figure 1.5: The Rotorcraft Noise Model (RNM) hemisphere measurement array. 


1.3.2.2 RNM/Q-SAM 

An extension to RNM enables it to model moderately accelerating or decelerat- 
ing quasi-steady flight conditions using the Quasi-Static Acoustic Mapping (Q-SAM) 
method, [13] which provides an aconstic equivalence between steady and quasi-steady 
straight-line flight conditions, allowing a relation to be made between the effect of 
moderate accelerations and changes in flight path angle. However, RNM is otherwise 
incapable of modeling the noise produced by unsteady flight and operating conditions 
other than those directly measured. The principle behind RNM/Q-SAM is discussed 


in more detail later, in Section 2. 2. 2. 2 
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Figure 1.6: The process of building a noise source hemisphere from 
measured data. 


1.3.2.3 HELENA and SELENE 


The European HELicopter Environmental Noise Analysis (HELENA) [H] and 
Sound Exposure Level starting from the Emitted Noise Evaluation fSELENEi|15j 
tools developed under the Friendcopter program are similar to RNM, except that 2D 
ground microphone arrays may be used. This allows noise to be measured across a 


wide range of directivity angles simultaneously. (Figure 1.8 on page 13 ) In concept, 


noise hemispheres could be developed from non-steady measurements, since the 
array covers nearly an entire hemisphere worth of directivity angles simultaneously. 
In practice, frequency spectra are generated at each microphone location from a 
hnite time window of acoustic pressure data. Depending on the required frequency 
bandwidth, this time window may be quite large — significantly larger than the time 
scale of acoustic events associated with the transient maneuver. Flight conditions may 


11 


Flight Path Angle, deg 


Airspeed, kts 

50 60 70 Noise Level 

^ ^ ^ > BVISPL 



Figure 1.7: RNM database of source noise hemispheres. 


be dehned with respect to either dimensional or non-dimensional parameters. For 
instance, dimensional parameters airspeed and rotor torque. [T3] and non-dimensional 
parameters advance ratio and tip path plane angle of attack, mi have been nsed in 
prior research. 

1.3.2.4 HAMSTER 

Qinetiq’s Helicopter Acoustic Measurement System for Trials and Experimental 
Reduction (HAMSTER) [3] technique is focused on developing models in order to 
assess the long-range detectability of helicopter operations. As snch, an accurate 
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Figure 1.8: The HELicopter Environmental Noise Analysis (HELENA) 
measurement array. 


description of the in-plane noise of the helicopter is desired. In addition to a ground- 
based array of microphones, HAMSTER uses tower-mounted beam forming arrays 
to measure noise in the plane of the horizon (Eigure 1.9 on the following page). 
By directing the array towards the helicopter, the influence of interfering acoustic 
reflections off of the ground can be suppressed. However, because of the narrow 
angle of the beam, the helicopter must be flown across the array at several altitudes 
in order to construct a full hemisphere. The relatively sparse linear microphone 
arrays have a limited number of microphone distance pairs; consequently, data is 
processed with a coarse 1/3-octave frequency resolution; constructing an accurate 
pressure time-history representation of the in-plane noise radiation would require a 
more rehned array. Due to the large amount of data required to construct a single 
model, hemispheres are typically only generated for several level flight speeds, and 
one climbing and one descending flight condition at a single moderate airspeed. 


1.9 on the following page 



Figure 1.9: The Helicopter Acoustic Measurement System for Trials and 
Experimental Reduction (HAMSTER) beamforming array. 
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1.3. 2. 5 Deficiencies of Source Noise Hemisphere Modeling 

While source noise hemisphere noise modeling is a major improvement over in- 
tegrated noise modeling, there are still significant deficiencies. The noise hemispheres 
developed under the compact source assumption cannot easily be generalized to other 
flight conditions, since each hemisphere represents the totality of noise radiating 
from the helicopter for a particular flight condition. While the hemispheres may be 
classified in terms of some non-dimensional parameters, the acoustic hemispheres 
still represent a single unique flight condition. As the flight condition changes, the 
contributions of each individual noise source to the overall external noise radiation 
changes. Consequently, a large flight test program is required in order to produce 
accurate estimates of helicopter noise radiation for a wide range of flight conditions. 
In general, the full range of helicopter operating conditions is not accounted for, 
especially variations in the configuration of the helicopter, such as changes in gross 
weight, and variations in the operating environment, such as air density or tempera- 
ture. Due to the expense of measuring the noise of the helicopter in all possible flight 
conditions, configurations, and atmospheric conditions, the effects of these variations 
are ignored in current helicopter noise source models. 
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1.4 Objective 


It is the objective of this dissertation to develop a new phenomenological 
method of modeling helicopter noise sources from measurements which: 

• Generalizes noise models based on the non-dimensional governing parameters 
of the noise sources 

In order to produce phenomenological models that describe the external noise 
radiation of the helicopter across a wide range of operating conditions, it 
is necessary to identify the governing parameters of the noise generation 
mechanisms. Multiple sources of noise contribute to the overall noise radiation 
of the helicopter; therefore, each source must be modeled individually with 
respect to its own governing parameters 

• Allows accurate noise estimates to be made at different operating conditions 
than those measured 

An accurate estimation scheme allows a wider range of operating conditions to 
be modeled with a smaller set of measured data. Operating conditions include 
not only various flight conditions, but also atmospheric and conhguration 
variations which affect the acoustic state of the vehicle. Current data-driven 
interpolation approaches are incapable of producing accurate noise estimates for 
operating conditions far from those measured. It is also impossible to accurately 
estimate noise radiation at operating conditions outside the measured range 
using interpolation techniques. Without the capability to estimate noise 
radiation at conditions not measured, general-purpose helicopter noise source 
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models will require an impractically large set of measured data to cover the 
full range of likely helicopter operating conditions. 

• Retains accuracy of noise estimates for measured flight conditions 

The new modeling method should preserve the accuracy of existing methods 
for flight conditions where measured data exists. The new method should 
categorically improve the quality of helicopter noise source modeling for ground 
noise contour estimation. 

1.5 Approach 

To achieve these objectives, a new modeling scheme is proposed. For the 
greatest flexibility and accuracy, each noise source on the helicopter must be modeled 
separately in terms of its own set of non-dimensional governing parameters. Existing 
physical knowledge of the major noise source mechanisms of the helicopter should 
be incorporated in the modeling structure, in order to improve the accuracy of noise 
estimates made for conditions not measured. The phenomenological models of the 
noise sources are then fit to specific flight vehicles using measured test data. The 
approach to developing this new modeling method is: 

1. Develop a method of separating the major noise sources of the helicopter. 

2. Develop a phenomenological modeling framework for each noise source in terms 
of non-dimensional governing parameters as well as a set of dependent modeling 
parameters which are used to capture physics not explicitly included in the 
assumed modeling structure 
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3. Employ parameter identification techniqnes to determine the valnes of the 
dependent parameters which cause the model to match measured data for flight 
conditions defined in terms of the non-dimensional governing parameters 

4. Validate the method using both wind tunnel and flight test measurements of 
rotor noise 


A diagram of the method is shown in Figure [T. 10 on the next page 
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Figure 1.10: Flowchart of the modeling approach developed in this dissertation. 
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Chapter 2 


Phenomenological Helicopter Noise Source Modeling 
2.1 Helicopter Noise Sources 

The comprehensive characterization of helicopter noise radiation is difficult. 
Helicopter noise is generated by many different sources. These sources can be broadly 
classified as rotor harmonic noise sources and non-rotor- harmonic noise sources. Rotor 
harmonic noise sources dominate the external noise radiation of helicopters. |16j Two 
types of rotor harmonic noise exist: lower harmonic noise sources and impulsive noise 
sources. The main and tail rotors continuously generate lower harmonic noise during 
all flight conditions. Impulsive noise only occurs during certain flight conditions, but 
is the dominant source of noise when it occurs. HU These noise sources are illustrated 
in Figure [2.1 on the following page 

Lower harmonic noise is produced by two mechanisms termed steady loading 
noise and thickness noise. Steady loading noise is generated by periodic application 
of aerodynamic forces by the rotor to the fluid medium — the resulting noise occurs 
at the blade passing frequency and its integer harmonics. For two-bladed and large 
multi-bladed main rotors, the first two harmonics will typically occur at sub-audible 
frequencies (i.e. below 20 Hz), observed only by electronic detection equipment or 
felt by induced structural vibrations. The higher harmonics of steady loading are of 
lesser magnitude and are generally significant only in the absence of more powerful 
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Figure 2.1: Helicopter rotor harmonic noise sources. 


noise sources. Thickness noise is generated by the displacement of the air caused 
by the periodic passage of the rotor blade sections through space. This effect is 
dominant near the tips of the blades where the section free-stream Mach number is 
greatest. The resultant noise is observed as a series of low-frequency pulses occurring 
at the blade passage frequency of the rotor. As the tip Mach number is increased, the 
intensity of thickness noise increases rapidly. The shape and magnitude of the pulse 
changes signihcantly as the advancing tip Mach number enters the transonic regime 
and local shocks develop about the blades — at this point it becomes a dominant 
source of noise and is termed High Speed Impulsive (HSI) harmonic noise. [18] HSI 
noise is dependent not only on the shape of the blade section but also on the transonic 
flow held surrounding the surface of the airfoil near the blade tip. HSI noise is less 
important on modern helicopter designs featuring lower tip speeds and thinner blade 
sections. 

The other impulsive noise source, termed Blade- Vortex Interaction (BVI) noise. 
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Frequency, Hz 


Figure 2.2: Typical frequency spectrum of helicopter noise. (Bell 206B3) 


occurs when the rotor blades pass near the vortices created and trailed from the tips 
of preceding blades. These vortices induce rapid changes in the local angle of attack 
of the following blades causing rapid changes in rotor airloads. Since the magnitude 
of loading noise is related to the time-rate-of-change of the forces exerted by the blade 
on to the medium, a highly impulsive noise results. This condition commonly occurs 
during moderate speed shallow descents. [IH] BVI noise is a particularly important 
consideration for community land use planning purposes because typical landing 
approach trajectories incorporate the same shallow descents which lead to high levels 
of BVI noise. Tail rotor harmonic noise sources are similar to main rotor harmonic 
noise sources and can be modeled in the same fashion. Typically, the tail rotor does 
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Figure 2.3: Notional acoustic pressure time history of helicopter noise. 


not operate in conditions where significant impulsive noise occurs — however, due to 
the higher rotational rate of the tail rotor, typically six to seven times higher than 
the main rotor, tail rotor lower harmonic noise sources occur at higher frequencies 
than main rotor lower harmonic noise. (Figure |2.2[ ) Since the human ear is more 
sensitive to these frequencies, tail rotor harmonic noise can be the dominant source 
of noise for annoyance. IZDl A notional diagram of the acoustic pressure time-history 


produced by rotor harmonic noise sources is shown in Figure 2.3 


In addition to the rotor harmonic noise sources, there are several noise sources 
that are not harmonic with the rotor blade passage frequency. An aerodynamic 
source of non-rotor harmonic noise is rotor broadband noise. The rotor produces 
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broadband noise over the mid- to high-frequency range. Such noise is caused by 
random aerodynamic interactions with the rotor blades, due to turbulence generated 
by boundary layer effects, vortex shedding, flow separation, and re-ingestion of the 
rotor wake, as well as atmospheric turbulence. |2T] This source of noise always exists 
on the rotor, but is usually only noticeable in operating regions where the other rotor 
noise sources are not dominant. While this noise source is generally not significant 
for helicopter rotors, it can be important for high thrust, low tip speed conhgurations, 
such as small unmanned vehicles or wind turbines. In addition to the aerodynamic 
noise sources, non-rotor harmonic noise is generated by the gearbox and engine of 
the helicopter. Engine noise has both tonal and broadband components, whereas 
gearbox noise is largely tonal. However, both sources of noise tend to be of much 
higher frequency than the other helicopter noise sources, and are much more readily 
absorbed by the atmosphere. As such, mechanical noise has little effect on ground 
noise contours except at very low altitudes where the distances between source and 
observer are small. This dissertation focuses on the modeling of rotor harmonic 
noise sources, since they are the dominant external noise sources of helicopters and 
consequently set the ground noise contours resulting from helicopter operations. 
Both main and tail rotor harmonic noise sources have the same mechanisms and 
can be modeled in the same fashion. One common way to model rotor harmonic 
noise sources is with the Ffowcs Williams - Hawkings (FW-H) equation, [22] which 
describes exactly how noise is generated by arbitrary surfaces in motion through a 
medium. 
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2.2 Harmonic Noise Source Modeling 
2.2.1 Ffowcs Williams - Hawkings Equation 


The FW-H equation (Equation |2.1[) can be divided into three terms: monopole, 


dipole, and quadrupole (Figure 2.4 on the following page). 
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The monopole term models thickness noise by considering the rotor blade as a 
set of monopole mass sources and sinks which describe how the blades displace the 
medium. The dipole term models the mechanisms of loading noise, including BVI, 
as a set of aerodynamic dipole sources on the surface of the blades that describe the 
forces the blades exert on the medium. The quadrupole term includes the effects of 
complex noise sources inside a fluid volume surrounding the rotor blades — this is how 
the effect of the transonic flow field that causes HSI noise is modeled. Since HSI noise 
is not a major consideration for the types of helicopters and operating conditions 
that are the focus of this dissertation, the quadrupole term will be neglected and 
only the on-surface monopole and dipole terms will be retained. 
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Figure 2.4: The aerodynamic origins of the three types of noise sources 
in the Ffowcs Williams - Hawkings equation. 
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2.2.2 Non-Dimensionalization 

2. 2. 2.1 Non-Dimensional Governing Parameters 


Previous research, both theoretical and experimental. [23] has shown that the 
external noise radiation of a particular rotor design in steady-state flight can be 
uniquely determined from a set of four non-dimensional governing parameters. In 
this dissertation, these parameters are non-dimensionalized as the rotor advance 
ratio (/i), hover tip Mach nnmber (Mh), thrust coefficient (Ct), and mean inflow 
ratio (A). The advance ratio and hover tip Mach nnmber control the motion of the 
blade sections through the medium, and hence the free-stream Mach nnmber seen 
by the blades over each rotor revolution. In addition, these parameters control the 
epicycloidal pattern of the wake (Fignre [2.5[ ), and hence the interaction geometry 
of individual blade vortex interactions. The interaction geometry controls how the 
acoustic effect of BVI accumulates in phase, and hence its magnitude and directivity. 


as shown in Fignre 2.6 on page 29 The thrnst coefficient and inflow ratio are 


important parameters in determining the magnitude and distribution of aerodynamic 
loads over the rotor disk. The inflow ratio inflnences the angle-of-attack of the 
rotating blade sections and hence lift coefficient. In addition, the inflow ratio controls 
the “miss-distance” between the rotor blades and the trailed wake dnring BVIs 
(Fignre [2.7[ ), and hence their strength. Lastly, the thrnst coefficient controls the 
rotor blade loading as well as the strength of the trailed tip vortices, and hence 
the strength of BVI. Sometimes, different sets of four independent non-dimensional 
governing parameters are selected, but they are equivalent - for example, advancing 
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tip Mach number {Mat) niay be substituted for hover tip Mach number, as the two 


are related to each other by advance ratio. (See Equation |2.9[) The selection of the set 


of parameters used in this dissertation will be shown analytically in Sections |2.2.2.3 
and 12.2.31 



Figure 2.5: The “top-view” geometry of the wake, controlled by /i and 
Mh- The potential BVI locations on the advancing and retreating sides 
are label. 



Figure 2.6: A diagram of the phase collection process of Blade- Vortex 
Interaction (BVI) noise. 
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Figure 2.7: The “side-view” geometry of the rotor wake, controlled by A. 
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2 . 2 . 2. 2 Generalization 


An important application of non-dimensionalization is that it allows noise 
models defined uniquely in terms of the governing non-dimensional parameters values 
to be generalized across a wide range of equivalent dimensionally-dehned flight 
conditions. The Quasi-Static Acoustic Mapping method (Q-SAM)[T9] is a prior 


application of this principle to helicopter noise, as mentioned in Section 1.3. 2. 2 


Q-SAM allows an acoustic equivalence to be made between flight path angle and 
low to moderate longitudinal accelerations. The equivalence is derived from the non- 
dimensional governing parameters — nnder the assumption that gross weight, rotor 
RPM, and airspeed are held constant, three of the four non-dimensional governing 
parameters can be held fixed: thrust coefficient, hover tip Mach number, and advance 
ratio. Then, by applying a simple longitudinal force balance and simple momentum 
theory, it is shown that rotor inflow ratio and hence BVI “miss-distance,” can be 
held constant for a constant rotor tip-path-plane angle of attack, which is itself 
defined by the constant helicopter drag-to-weight ratio and the variable helicopter 


flight path angle and longitudinal acceleration. (See Section 2.2.3. 1 for more details.) 
Consequently, for moderate accelerations where the airspeed of the helicopter changes 
slowly, accelerations can be modeled by using measured data for steady-state flight 
conditions at a flight path angle resnlting in an eqnivalent tip-path-plane angle of 
attack, and hence eqnivalent inflow ratio. When combined with RNM, the effect is 
the selection of acoustic hemisphere models based on main rotor advance ratio and 


inflow ratio, as opposed to airspeed and rate of descent, as shown in Figure 2.8 on 
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Figure 2.8: RNM/Q-SAM database of source noise hemispheres based on /r and A. 


This concept was further expanded to steady-turning flight conditions by the 
author, [21] allowing moderate accelerating, turning, and descending maneuvers to 
be modeled using steady straight-line flight data. However, under steady turning 
flight conditions, thrust must vary to counteract centripetal acceleration. Variation 
in thrust changes the inflow ratio, the amount of which can be determined from 
momentum theory, but also changes the vortex strength for which the effect on noise 
cannot be directly determined without the acquisition of helicopter noise data at 
several thrust levels. The result is a helicopter noise model defined in terms of main 
rotor thrust coefficient in addition to advance ratio and inflow ratio, as shown in 


Figure 2.9 on the next page 
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Figure 2.9: RNM/Q-SAM database of source noise hemispheres, extended 
to steady turning flight. 


2. 2. 2. 3 Non-Dimensional Ffowcs Williams - Hawkings 

In order to facilitate comparisons between acoustically equivalent rotor operat- 
ing conditions, it is benehcial to rearrange the FW-H equation in a non-dimensional 
form. This can be performed, based on the approach used in Reference [23], as 
follows: 

First, the time of observation of the acoustic pressure is non-dimensionalized 
by the rotational rate of the rotor: 
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The time of emission is likewise non-dimensionalized: 


r 

n 


The geometric parameters can be non-dimensionalized — following standard 


convention, these are non-dimensionalized by rotor radius: 


T 
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(2.4a) 


(2.4b) 


(2.4c) 


The Mach number is dehned as the ratio between the speed of the element 
with respect to the medium to the speed of sound: 


ttQ 


The hover tip Mach number is defined as the relation between the rotor radius 
and the rotational rate: 


Mm = 


The tip Mach number in forward flight can then be dehned in terms of the 
advance ratio and related to hover tip Mach number: 


U QR + Vsini; 

Mt = — = = Mht (1 + /ismr) 

Oq «o 


where the advance ratio is the ratio between the aircraft velocity and hover tip 


velocity: 
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The advancing tip and hover tip Mach numbers are therefore related by advance 

ratio: 

Mat = Mht (1 + (2.9) 

By neglecting spanwise and perpendicular flow, the velocities normal to the 
rotor blade surfaces are assumed to have the form: 

dz 

Vn = U— = Ui = Maoi (2.10) 

dy 

where ^ is the geometric slope of the airfoil surface. 

The acoustic pressure is then non-dimensionalized as the acoustic pressure 
coefficient: 

= (2.11) 

Po% 

Similarly, the on-surface blade pressures may be non-dimensionalized with 
respect to the dynamic pressure: 


- 1^0 


(2.12) 


Substituting Equations 2.2 through 2.12 back into the FW-H equation (2.1), 


and neglecting the quadrupole term yields: 
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Simplifying this expression yields a non-dimensionalized version of the FW-H 


equation (2.1): 
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Neglecting the spanwise and perpendicular flow components, and substituting 


for the resulting blade element Mach numbers: 
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where 6 is the angle between the assumed rotor plane and the observer. 


From Equation 2.15 it is apparent that monopole thickness noise can be solved 


for directly with operating condition of the rotor defined by hover tip Mach number 


(Mh), advance ratio (/i), and the blade geometry. For higher tip Mach numbers, 
this formulation is known to underpredict the resulting noise levels due to additional 


off-surface aerodynamic effects which are included in the neglected quadrupole 


volumetric source term. However, for moderate tip Mach numbers, the monopole 
prediction provides a reasonable level of accuracy. 

The dipole term is likewise dependent on hover tip Mach number (Mh) and 
advance ratio (/i); however, it is also a function of the non-dimensionalized on-surface 
blade pressures described by Cp^.. Further assumptions about the aerodynamic 
loads must be introduced to identify the complete set of non-dimensional governing 
parameters for rotor harmonic noise. 


2.2.3 Aerodynamic Loading Model 

Solving for the dipole loading noise of the rotor is a more complex task than 
the monopole thickness noise calculation. While the monopole source strengths are 
known exactly from the geometry and motion of the rotor blades with respect to the 
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medium, the determination of the dipole source strengths is less straighforward. The 
four non-dimensional governing parameters of rotor harmonic noise (/i, Ct, A, and 
Mh) specify a unique operating condition for the rotor, but there is no simple way 
of calculating the exact distribution of loads across the rotor in both space and time. 
Without restrictions on how loads can be distributed across the rotor disk, there are 
a theoretically inhnite number of conhgurations that can yield the same far-held 
noise radiation. However, by introducing some reasonable engineering assumptions, a 
physical structure can be developed to limit the set of possible loading distributions 
to those that are practically realizable. 


2.2.3. 1 Trim 


The hrst component to the assumed loading structure is a simple trim model. 
Given an estimate of the hat plate drag area and rotor “hub-force” (related to prohle 


drag) of the helicopter, a longitudinal force balance (Figure 2.10 on the following 


page ) can be established in order to estimate the angle of attack of the tip-path-plane 


for a given hight condition as follows: 


Df + H ttx 


where the parasite drag-to-weight ratio can be estimated as: 

Df _ fe 


(2.16) 


W 2ACt 


e 2 
/i 


(2.17) 


For low hight speeds, H/W may be included in the parasite drag term, otherwise 
it may be estimated using a simple blade-element theory derived approximation, one 
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standard expression being: 


(2.18) 


H acdo 1 + 4.6/i^ 

W 8 



Tw 

Figure 2.10: Longitudinal force balance. 


From the angle of attack of the rotor disk and the required rotor thrust, the 
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non-dimensional inflow throngh the rotor can be calcnlated numerically using the 
equation: 


A 


— fjj tan O.TPP ~\~ 


Cp 

2yA2 + p 


(2.19) 


where the inflow ratio is deflned in forward flight as: 


A 


—V sin appp + Vi 
OR 


( 2 . 20 ) 


and the thrust is non-dimensionalized as thrust coefficient, which is defined as: 


Cp — 


T 

PoAinny 


( 2 . 21 ) 


2. 2. 3. 2 Lower Harmonic Loading 


Next, the lower harmonic loading can be given an assumed form. Under the 
assumption of uniform inflow and an ideally twisted rotor, blade element theory 
suggests that the spanwise distribution of rotor blade loads should be approximately 
triangular from the blade hub to the blade tip. In the chordwise direction, the load 
is concentrated at the quarter chord. The azimuthal distribution can be specified in 
rotor time for each individual lower harmonic using a truncated Fourier series: 

k 

T = To + Tnc cos ml; -\- Tns sin nip (2.22a) 


D = Do 


Dnc COS nip -\- Dns sin nip 


(2.22b) 


n 


A diagram of the T\c contribution to the loading distribution is shown in 


Figure 2.11 on the next page 
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Figure 2.11: Contribution of Fourier coefficient T\c to lower harmonic loading. 


2 . 2 . 3. 3 Rotor Wake Modeling 

While suitable for the lower harmonic loads, a direct specihcation of the 
individual loading harmonics is tedious for impulsive noise sources, such as BVI, 
and does not include additional physical understanding of the problem. A better 
approach is to model the BVI noise process with a model of the vortex wake structure, 
providing a physical description of the blade-vortex interaction, including both the 
effects of the miss-distance between the vortex and the rotor blade, which determines 
the strength of individual interactions, as well as the phase collection process of the 
interactions, which determines the directivity of radiated noise. 


Figure 2.12 on page 42 uses Hugens’ wavelets to illustrate in two-dimensions 


how the interaction angle between the vortex and the blade causes the acoustic effect 
of the BVI to phase in the medium causing the noise to be radiated in different 
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directions. For instance, on the advancing side of the rotor the “a” interaction 
occurs further from the rear of the rotor and at an oblique angle; as the disturbance 
moves along the interaction through the medium, the waves collect in phase focusing 
the radiated noise ahead of and towards the retreating side of the rotor. As BVI 
occur farther back on the rotor disk, the interaction angle changes and directivity 
shifts towards the advancing side. Interaction “ 7 ” is an example of a parallel BVI, 
where the blade and vortex are parallel at the time of the interaction so that the 
disturbance in the medium occurs near-simultaneously along its extent-the results 
in a large magnitude of BVI noise radiated towards the advancing side of the rotor. 
The same process occurs on the retreating side of the rotor, along these interactions 
contribute less to the overall noise radiation of the helicopter since the blade section 
Mach numbers are lower on this side. 

The BVI noise process is dominated by the trailed rotor tip vortices; in order 
to simulate BVI, a physics-based description of the path of the trailed tip vortices 
through medium is desired. Undistorted rigid wake models lack the flexibility 
required to accurately model the miss-distance of the interactions. On the other 
hand, distorted prescribed wake models provide a realistic structure for describing the 
path of the trailed tip vortices through the medium for a given flight condition. There 
are two major prescribed wake models in use: the UTRC generalized wake mo del [25] 
and Beddoes wake model. [26] The generalized wake model defines a set of geometric 
distortions to a rigid wake. In practice, the values of these distortion parameters 
are determined from empirically measured data for specihc flight conditions and 
rotor conhgurations. The distortions to the wake in this model have no direct 
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Figure 2.12: Acoustic phasing and corresponding in-plane directivity of 
different advancing and retreating side BVI. [1] 


physical meaning and are simply a description of the measured wake geometry. The 
Beddoes wake model, on the other hand, dehnes tip vortex locations for steady 
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flight under some physical assumptions. In the original Beddoes wake model, the 
position of the tip vortices was determined by assuming a linear longitudinal inflow 
distribution superimposed with a lateral cubic downwash distribution to account for 
tip effects. The convection of the non-interacting vortices is then calculated based 
on the combined freestream flow and prescribed induced inflow distribution over the 
rotor. That is: 

ri’v 

z = — Xid'ip (2.23) 

J'<Pb 

where 0 = V’b ~ f/’u is the wake age and where the Glauert[2^ longitudinal 

linear inflow distribution is combined with a lateral cubic downwash distribution to 

produce Beddoes’ prescribed inflow over the rotor disk such that: 

\i = Xig {1 + X cos ij + x\rv sin'll; f) (2.24) 


The mean induced inflow ratio, can be solved for numerically from 


Aio 


Cj 


‘^\j if^z + Aj)^ -I- 
X is the wake skew angle, defined as: 


(2.25) 


X = 


tan 


-1 


/^z + A 


*0 


(2.26) 


The Beddoes wake method can be further extended by assuming alternate 

inflow distributions across the rotor disk. For instance, van der Wall [29 j incorporated 

Drees’ [30] lateral inflow model and increased the constant inflow ratio over the area 

swept by the rotor blades in order to account for the effect of the root cut-out: 

Xi = Aig (1.2 -f ycosV’ + x| r^sin^’l^ - 2/ia,sin^/>) (2.27) 


In practice, the best coefficients of the linear inflow distribution are not known 
a priori for a given rotor configuration. For instance, the longitudinal inflow ratio 


43 


is sometimes taken as | instead of x, based on other experimental observations of 
rotor inflow. A more general way of expressing a range of possible assnmed inflow 
distribntions is: 

K = Ko {^X (cos -0+1 sin -01^) + + C) (2.28) 


where A, B, and C are parameters which control the respective longitndinal, 
lateral, and mean components of the prescribed linear distribution of inflow over the 
rotor. 


Given this general form of the assumed inflow distribution, the vertical displace- 


ment of the rotor tip vortices can be solved for using Equation 2.23 Underneath 


the rotor disk, where BVI occur, the integration results in the following algebraic 
expression of the vertical wake geometry: 


z = — A, 


*0 


Ax ^cos-^^ -f — \vy sin'^^l^^ -|- 2Bjjix sin'^^ -|- C 


0 (2.29) 



Figure 2.13: A side view of the Beddoes wake model. 


In addition to the vertical displacement of the wake, the longitudinal and 

lateral motion of the wake affects the phase collection process of the BVI. Near the 

rotor disk, the “top-view” geometry of the wake is mostly epicyclodial, although the 

tip vortices roll-up somewhat inboard of the blade tips: 

X = cos'ipv + (2.30a) 
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Lateral Wake Skew 



Figure 2.14: A head-on view of the Beddoes wake model. 


y = r^ sin 


(2.30b) 



Figure 2.15: An overhead view of the Beddoes wake model. 


Over time the rotor wake contracts slightly, which can change the phase of late 
BVI relative to the earlier BVI. This can be accounted for by incorporating a model 
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(2.31a) 

(2.31b) 


of the wake contraction, such as that identihed by Landgrebe. [31] 

ry = Tr [D + {1 - 

Ai = 0.145 + 27 Ct 


In addition to the wake geometry, the structure of the vortex has an effect on 
the magnitude and impulsiveness of the interaction. The Biot-Savart law describes 
the flow induced by an inviscid vortex, and accurately models the flow induced by 
real vortices at relatively far distances from the vortex. However, near the vortex, 
viscosity reduces the velocity induced by the vortex. The Vatistas model[32| describes 
the resulting tangential velocity prohle induced by an inhnitesimal viscous vortex 
segment: 


_ IvT 

^ 27r(r^”' T 7 ’ 2 n^(l/n) 

If the vortex strength is non-dimensionalized as: 


(2.32) 


Iv = 


Iv 


(2.33) 


Then the model can be expressed as: 


ve IvT 


VLR 27r(r2"' + 7-2nyi/n) 


(2.34) 


where r = r/R and = rjR. 

The Vatistas model is a generalization of several commonly used vortex models 


(Figure 2.16 on the next page). For n = 1, this is equivalent to the Scully vortex 
model. [33] For n = 2, the model is equivalent to the Bagai-Leishman model. [M] In 


the limit n — >■ cx), the Rankine vortex is described. In this dissertation, the commonly 
used Bagai-Leishman form with n = 2 is assumed, although the choice of n has little 
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impact on the resulting BVI noise since the slope of the tangential velocity inside 
the core region is similar for all values of n. The viscous vortex core does not stay a 
constant size, but grows as the vortex ages. One common analytical model of viscous 
vortex core growth is the Lamb-Oseen vortex model. [35] In non-dimensional form, 
this can be expressed as: 

f, = y/4^ (2.35) 


where is an empircal constant. 



Figure 2.16: Vortex core tangential velocity prohles for several common 
models, normalized against Rankine core velocity. 


However, the Lamb-Oseen model starts with an inhnitesimal vortex core size. 
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which is not physically realistic. Bhagwat and Leishman[36] modihed the Lamb- 
Oseen model to inclnde an initial vortex core size, which more closely matches 
experimental measurements. 

fc = \jrl + ACy(f) (2.36) 

The non-dimensional strength of the vortex is assumed to vary harmonically 
about the azimuth: 


r = To (70 + 7i5 sin + 7ic cos -0^) 


(2.37) 


where the nominal value of the vortex strength is that produced by constant spanwise 
linear blade loading, again non-dimensionalized by flR^: 


Bn = 


2tiCi 


(2.38) 


Once the wake structure is fully defined, the instantaneous velocities induced 
on the blades as they travel about the rotor disk can be determined by summing the 
contributions of the trailed tip vortex segments, and from these the time-varying 
angles of attack of the blade sections. 

Due to the relative simplicity of this wake model, it may not be suitable for all 
rotor configurations. The Beddoes wake model is a good choice for two-bladed rotors, 
but may not adequately describe the more complex wake geometry of multi-bladed 
rotor systems. Likewise, the single trailed tip vortex assumption is acceptable for 
lightly twisted rotors, but is known to be inadequate for highly twisted or flapped 
rotors. The modeling methodology developed in this dissertation will be applied 
to linearly-twisted two-bladed rotors in Chapter where the assumed wake model 
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is known to be valid. However, the assnmed wake model may be inadeqnate to 
accnrately characterize BVl for other rotor conhgnrations and mnst be validated 
and rehned for these conhgnrations, if necessary. 


2. 2. 3. 4 Unsteady Aerodynamics 

In the assnmed wake model described in the previons section, the presence 
of the shed wake on the trailed wake geometry is assnmed small and has been 
neglected. However, dne to the shed wake, the effect of the nnsteady velocity 
indnced on the rotor blades by the tip vortices is not felt instantaneonsly and this 
effect shonld be incorporated in the BVl model. A common way to acconnt for 
nnsteady aerodynamic effects is throngh an indicial aerodynamics model. As this 
dissertation considers only hight conditions at low to moderate advancing tip Mach 
numbers, the Leishman-Beddoes 2D incompressible indicial aerodynamics model has 
been chosen. The Duhamel integral is evaluated over time at each blade station to 
determine the effective compact chord aerodynamic angle of attack, using the mid- 


point integration algorithm developed by Beddoes|3H](Eq. 2.40) with the empirical 


constants (Ai,A 2 ,&i,& 2 ) determined by Leishman. 

Q!e('S) = a(s) — A(s) — V (s) 


(2.39) 


where. 


A(s) = X{s - 

V(s) = V(s - As)e-^2^" + 


(2.40a) 

(2.40b) 


where s is the reduced time, i.e. the distance traveled by the airfoil in terms of 
semi-chords over a given time interval t: 


s_2 


Vdt 


( 2 . 41 ) 
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Given the unsteady effective angle of attack for each blade section, the section 


lift and drag can be readily calculated from tabulated airfoil section data, and from 
those values the non-dimensionalized surface pressures, Cp-y The acoustic pressure 
time-history can be readily calculated for any observer using Farassat Formulation 


f A|40j of equation 2.f4, thereby providing a complete analytical model of rotor 


harmonic noise generation. 


2.2.4 Choice of Dependent Parameters 

Altogether, the assumed loading model contains many empirical parameters 
that vary depending on the flight condition of the rotor as determined by the non- 
dimensional governing parameters of loading noise as well as the specihc conhguration 
of the rotor. Given a set of values for the unknown dependent parameters and the 
known independent governing parameters, the noise radiation of the helicopter is 
uniquely determined. However, the association between the dependent parameters 
and the governing parameters is not known a priori. The values of the dependent 
parameters must therefore be identified for each set of non-dimensional governing 
parameters based on measured data for a particular rotor conhguration. Altogether, 
the set of dependent parameters to be identified is: 

X = ^ A B C D Tr ro Cv JO 7ic 7is To Tic Tis ■■■ Do Die Dis ■■■ ]■ (2.42) 
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Chapter 3 


Parameter Identification 

Parameter identification is a method of building a model by determining values 
for a set of parameters, which are unknown a priori, in order to match the output 
of the model to measured data. Determination of the unknown parameter values 
requires solving an inverse problem, i.e. given the desired output for the model, 
determine a set of input parameters to the model. In many cases, analytical solutions 
to the inverse problem are not available and numerical solutions must be pursued. 
Numerical solutions to inverse problems iteratively run the assumed model in the 
forward direction for various trial parameter values until a combination is found 
which results in the model output closely matching the measured data. To achieve 
feasible solution times, this requires that the assumed model can be evaluated quickly. 
In addition, the assumed model should have no more degrees of freedom than are 
required to match the measured data — this ensures that the space of trial parameter 
values is not too large and also reduces the likelihood that there will be multiple 
distinct and incompatible solutions to the problem. 

The parameter identihcation problem can be divided into three components. 
The hrst is the assumed analytical model of the phenomena to be modeled. The 
second is a quantitative measure of how closely the model output matches that 
measured data, termed an error metric. Lastly, an efficient method of choosing 
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trial parameter values so as to reduce the error metric, an optimizer, must be 


chosen. Altogether, these components set up the parameter identihcation method 


as an optimization problem, where the optimal values of the unknown dependent 


parameters are determined so as to minimize the error between the model output 
and measured data. A flowchart diagramming the complete method is shown in 



Figure 3.1: Flowchart of the parameter identihcation process. 
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The physics based modeling structure for rotor harmonic noise dehned in the 
previous section provides the basis for the assumed model used in the parameter 
identihcation technique. The unknown parameter values to be identihed are the 
dependent parameter values from the assumed loading model, for a given set of 
known non-dimensional governing parameters. The ontpnts of the model are the 
acoustic pressure time-histories produced over one rotor period at a set of observer 
locations surrounding and traveling through the medium with the rotor. 

3.1 Error Metrics 

3.1.1 Time-Domain Error Metric 

The error metric describes how closely the model output matches measured 
data. In conjnnction with the assnmed physical model, this forms the objective 
function for the optimization problem posed by parameter identihcation. The error 
metric conveys information abont the qnality of the trial parameter valnes to the 
optimizer. The best choice of error metric does not necessarily have the same form 
as the metric that will be used to quantify the output of the identihed model. The 
rotor acoustic models developed in this dissertation will be used to estimate ground 
noise contours, typically quantihed using a frequency-weighted Sound Pressure Level 
(SPL); for example, the A- weighting is often nsed to model hu m an perceptions 
and annoyance. However, these metrics rehect only the total energy content of the 
signal — there can be many acoustic signals with the same SPL but radically different 
waveforms, and hence physical mechanisms. In order to develop models which more 
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closely match the actual mechanisms of rotor harmonic noise, the error metric shonld 
include contributions of both the frequency and phase of the signal. This can be 
accomplished by formulating the error metric in the time domain. One common way 
of doing this is to use the generalized ;?-norm, which is an integrated measure of error 


over an interval of time, in this case one rotor blade passage, as shown in Figure 3.2 

1 

rT 


\e\\p = 




(3.1) 


The value of p describes how the greatest error at any point is weighted against 
the average error over the entire time period. A value of p equal to unity is equivalent 
to the average error over the signal. A value of p equal to two is the RMS error of 
the signal. In the limit, as the value of p — )■ oo, the norm equals the maximum error 
at any point in the interval. 

The generalized p-norm describes how the error metric may be formulated 
for a single observer, however the pressure time history at one observer location is 
not unique and could be achieved by several distinct sets of parameter values. The 
uniqueness of the parameter identihcation solution can be improved by including 
additional information about the problem in the chosen error metric. Rotor harmonic 
noise is highly directional, so including information for several distinct observer 
directions in the error metric will aid in determining a single physically realistic set 


of dependent parameter values. (Fig. 3.3) 


54 


Time Domain Observer Error Metric 



Figure 3.2; Calculation of the time domain integrated error metric. 
The solid black curve represents measured data, the blue dashed curve 
simulated data. The error metric is an integrated measure of the absolute 
value of the difference between the two. 


3.1.2 Multi-Observer Error Metrics 


An obvious way to formulate an error metric for several observer locations 
is to take a weighted sum of the individual time-domain error metrics for each 


observer location. (Eq. 3.2) For a perfect measurement and model, the choice of any 


non-zero weighting does not influence the solution, since all individual error metrics 
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can be simultaneously reduced to zero. However, with the addition of measurement 
noise and modeling error, it will be impossible to achieve zero error at all observers 
and therefore a Pareto-frontier of solutions will exist in the parameter space. The 
Pareto-frontier is the surface of solutions where the error cannot be reduced at any 
one observer without increasing it at one or more others. In this case, the choice of 
observer weightings will determine which directions will be more accurately modeled, 
and which less accurately modeled. The choice of weighting can be tailored to the 
intended use of the noise model — for instance, a model for long range detection 
may weight the in-plane observer locations more heavily than out-of-plane observer 
locations, vice versa for a community noise model. 

^obs 

F{X) = J2^k\\e\\’; ( 3 . 2 ) 

k=l 
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Another approach is to formulate the problem as a multi-objective optimization 
problem. In this case, instead of attempting to find the best solution for a single set 
of weights determined a priori, the optimizer attempts to define the entire Pareto- 
frontier, allowing the determination of observer weighting to occur after analysis. 
This multi-objective optimization approach requires a more expansive search through 
the parameter space, and hence a higher computational cost. 

In addition to relative weightings of the observer locations, the number and dis- 
tribution of observer locations to be included in the error metric must be determined. 
A sufficient number of observer locations must be chosen in distinct directions so as to 
discriminate the different rotor noise sources based on their directivity characteristics. 
However the availability and uncertainty of measurements in certain directions will 
also dictate the choice of observer directions included in the metric. In addition, the 
computational cost of evaluating the model will increase linearly with the number of 
observers 

3.2 Optimization Methods 

3.2.1 Conventional Optimization 

Once the model and error metric have been determined and combined to 
formulate an objective function, the parameter values that minimize the objective 
function must be determined. An enumeration of all possible parameter values 
throughout the solution space would eventually find the optimal solution, but the 
computational cost of doing so would be inordinately high. A better approach is 
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to use optimization techniques to make an informed choice for each new set of trial 
parameters based on previous results. Traditionally, gradient-based optimization 
techniques are employed to minimize the objective function. These techniques start 
at an initial guess for the parameter values. At this point, the gradient of the 
objective function is calculated, indicating the direction from the trial point in which 
the error is reduced. If the gradient is not known analytically, the objective function 
will have to be evaluated at several points in the vicinity of the trial point in order to 
calculate the gradient by hnite difference. The objective function is then evaluated 
several times in a direction determined by the gradient, so as to lead to a reduced 
value for the objective function. The set of parameter values dehning this point 
in the parameter space become the new trial point. The process is repeated until 


no further progress is made in reducing the objective function value. Figure 3.4 
shows the objective function (i.e. error metric) contours for a simulated set of lower 
harmonic data — the third harmonic thrust dimensions are shown as representative 
examples. This type of objective function is smooth and convex, ideally suited to 
conventional optimization techniques. 


Figure ^ shows the simulated objective function contours for the model with 
both lower harmonic and BVI noise. These contours are shown in terms of two 
parameters which control the wake: the longitudinal wake skew coefficient (A) and 
the wake contraction coefficient (D). Unlike the lower harmonic only objective 


function shown in Figure 3T the inclusion of the wake model leads to a non-convex 
objective function with multiple local minima. In cases where the objective function 
contains numerous local minima, gradient-based methods will halt at the hrst local 
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2D Objective Function Contours 



Figure 3.4: Simulated third harmonic thrust coefficient lower harmonic- 
only objective function contours. 


minimum encountered, which may not be similar to the global optimnm over the 
entire parameter space of the problem. 

3.2.2 Global Optimization 

In order to determine the global optimum of a non-convex objective function, 
a global optimization scheme must be employed. Stochastic methods are designed 
to explore the entire parameter space and identify probable global optima with a 
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2D Objective Function Contours 



Figure 3.5: Simulated wake parameter objective fu n ction contours for 
lower harmonic and blade-vortex interaction noise model. 


minimum number of function evaluations. Since these methods explore the entire 
parameter space, they are also well suited to identifying the Pareto-frontier in multi- 
objective optimization problems. Stochastic methods employ the statistical properties 
of random number distributions to ensure that the solution space is explored and 
that the method does not terminate in potential local minima. Common stochastic 
methods include genetic algorithms, HU differential evolution. |42j and more recently, 
particle swarm methods. |43j Genetic algorithms are well suited to non-convex and 
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non-smooth problems, and are capable of identifying global optimnm over a wide 
and varied parameter space. The genetic algorithm starts by seeding the parameter 
space with a population of points at randomized locations. The parameter values of 
each member of the population are then discretized and encoded into a binary form. 
The binary values of each parameter are concatenated into a single binary number 
representing the particle location, called a chromosome. A subset of the population is 
paired with other members, and a randomized portion of the chromosomes swapped 
to form new chromosomes. Then, random bits of the chromosome are mutated by 
switching bits from 0 to 1 or vice versa. The chromosomes are decoded back to 
parameter values, and the objective function evaluated for each. Then, the members 
of the population are removed through a roulette selection, with the probability 
of a member being removed increasing as in proportion with the associated error 
value of the objective function. Over successive iterations, the members of the 
population with the lowest error values tend to survive and share information with 
other members of the population. However, the chromosomes of the members of the 
genetic algorithm relate only indirectly to the contours of the objective function over 
the parameter space. This lack of a direct relation between the parameter space and 
the chromosome values allows the genetic algorithm to explore the entirety of the 
solution space in a robust manner, even in the presence of discontinuities. However, 
for many physical optimization problems, the objective function may be smooth 
over the entirety of the parameter space, and while multiple local minima may exist, 
they are likely confined to a small region. On such problems, the genetic algorithm 
spends many function evaluations exploring in fruitless directions. 
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The differential evolution algorithm is an adaptation of the genetic algorithm 
that defines each member of the population more directly by the vector of its 
parameter values. Parameter values are swapped between members of the population 
during pairing, and during mutation, the parameter values are perturbed. In this 
way, differential evolution can more directly exploit the shape of the objective 
function across the solution space, and search for optima where they seem more 
likely to occur. However, the amount of perturbation is governed by several user 
adjustable parameters. The ability of the differential evolution method to find the 
global optimum is very dependent on the choice of these parameters, which can be 
substantially different for different objective functions. Without careful tuning of 
these parameters, the likelihood of missing the global optimum altogether is high. 

A more recent approach which has been found effective for physically-based 
non-convex objective functions is called Particle Swarm Optimization (PSO).[13] 
PSO evaluates the objective function at a number of randomly seeded particle 
locations in the feasible solution space. The particles are assigned velocities dictating 
how they travel through the solution space. At hxed intervals, the objective function 
is again evaluated at each particle location, and a new velocity computed for each 
particle. The velocity of a particle during the next iteration is a randomly 

weighted combination of the current velocity direction (vj), the direction towards 
the previous best set of parameters for the particle and the global best set of 

parameters discovered by any particle The common PSO update algorithm 

can be expressed simply: 

!.'+■ = av‘ + (.^1 - X‘) + hfl), (X‘, - X‘) (3.3a) 


62 


Xl + vl 


(3.3b) 


Xf ^ = 

where and ^2 are random variables between 0 and 1, a is the inertia weight of 
the particle, and bi is the weight of the direction towards the historical best position 
of particle i, and 62 the weight of the direction towards the best position of the entire 
swarm. There are several common sets of weighting terms for the PSO algorithm 
that are broadly applicable to a wide range of problems — in this dissertation, the 
weights suggested by Trelea are used. |H] 



Figure 3.6: The Particle Swarm Optimization (PSO) technique. 
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The stochastic nature of the algorithm encourages the particles to explore 
the entire parameter space while the inertia term allows the particles to pass over 
local minima. However, the quick convergence of the standard PSO approach still 
makes it susceptible to premature stopping. Several extensions of the method have 
been proposed to improve the ability of PSO to hnd the global optimum. One 
approach employed in the PSO method used in this dissertation is the “particle reset” 
method. |in] When the velocity of a particle decreases below a small predetermined 
threshold, the particle’s location is reinitialized to a random position in the space 
of possible solutions. This method prevents individual particles from stagnating in 
local minima, and ensures that the method continues to spend some time searching 
throughout the entire search space, even after a promising minimum has been found 
by the swarm. 

All three global optimization methods were applied to the parameter iden- 
tihcation problem for rotor harmonic noise using simulated data. In general, the 
PSO optimization consistently produced good results without much hand tuning of 
the optimization parameters after a relatively small number of iterations, and was 
selected as the optimizer for the method developed in this dissertation. However, it 
is likely that good performance can be achieved with any of the global optimization 
methods once they are adapted to this specific parameter identification problem. 

Another method to improve the ability of global optimization methods, includ- 
ing PSO, to hnd the global minimum is to initialize the method so as to provide 
better coverage over the entire search space. The simplest way to improve the initial 
coverage of the search space is to increase the number of particles in the swarm, but 


64 


this increases the computational cost of the method in proportion to the size of the 
swarm. However, coverage can be improved without impacting the computational 
cost. The standard approach of seeding the method with a collection of particle 
locations determined using a uniform random number generator results in a very 
non-uniform distribution of particles over the solution space, i.e. the seeded particles 
tend to cluster in certain areas of the solution space and are sparsely populated 
in others. The uniformity can be improved by using a “quasi-random” number se- 
quence, which defines streams of numbers which are more uniformly spread through 
the solution space, and yet still possess many of the statistical properties of truly 
random numbers. Numerous quasi-random number sequences exist — one which has 
been successfully applied to Monte Carlo methods [Mj and stochastic optimization 
algorithms [TF]. including PSO[l8j, is the Sobol sequence. [19] The Sobol sequence is 
used in this dissertation to seed an initial population of particles numbering ten times 
the dimensions of the solution space. This “quasi-random” seeding allows a smaller 
number of particles to ensure adequate coverage of the solution space, leading to a 
higher likelihood that the global minimum will be correctly identihed. In this disser- 
tation, 200 “quasi-randomly” seeded particles are used to provide complete coverage 
over the dependent parameter space during the PSO parameter identihcation. 

3.3 Parameter Estimation 

The parameter identihcation technique provides the means to develop an 
association between a set of governing parameter values and dependent parameter 


65 


values for a single operating condition. To estimate noise for an operating condition 
not measured, it is necessary to determine the set of dependent parameter valnes 
that correspond to the known governing parameter values which dehne the operating 
condition. 

3.3.1 Conventional Interpolation 

One method is to use an interpolation method to estimate dependent parameter 
valnes based on dependent parameter valnes identihed for other sets of governing 
parameter values. There are several practical interpolation approaches. The simplest 
method is the nearest neighbor approach, where the dependent parameter valnes 
are simply chosen to match the closest available measnrement. This method is 
suitable when known data points are tightly spaced, and estimates only required at 
positions similar to measured data. Linear interpolation allows for estimates of the 
dependent parameter valnes to be made for sets of governing parameters that are 
farther removed from known values. There are two possible methods of performing 
linear interpolation on the mnlti-dimensional space of governing parameters. The 
simplest is to perform 1-D linear interpolation separately along each dimension in a 
predetermined order. This is similar to the approach currently utilized by RNM to 
interpolate the data values on aconstic hemispheres as a function of both airspeed 
and flight path angle. This approach ties the accnracy of the interpolation to the 
choice of governing parameters and the order in which the ID interpolations occur. A 
better approach is to perform an iV-dimensional linear interpolation, where N is the 
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number of governing parameter dimensions. To accomplish this process, the space of 
governing parameter values must be triangulated using the measured parameter sets 
as vertices. Then, dependent parameter values can be linearly interpolated along the 
hyperplane formed by the nearest set of measured value vertices. This method still 
has some geometrical accuracy dependence, requiring data be gathered uniformly 
across the governing parameter set in order to produce good results. In addition, 
this method cannot yield extrapolated values outside the range of measured data. In 
order to make the best use of measured data, a more sophisticated approximation 
scheme capable of extrapolation and able to yield better results with scattered data 
is sought. 

3.3.2 Radial Basis Function Approximation 

One increasingly common way to approximate unknown nonlinear functions 
is the Radial Basis Function (RBF) approach. |5Uj Any smooth function may be 
approximated as a linear combination of radial basis functions. Each RBF is 
centered at some distinct point in space. Each RBF is weighted — these weights are 
selected to achieve the best-known approximation of the target function. The RBF 
approximation can be expressed in the following form: 

N 

f{x) ^ ^ Ci)) (3.4) 

i=l 

The radial basis functions, $, are weighted so as to match the dependent value 
at each set of independent coordinates. The argument, r, of the radial basis function 
is some scalar representation of the distance between the approximated point, x, and 
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the center of the radial basis function, q. In rectilinear coordinates, this distance is 
generally taken to be the Euclidean distance between the two points, and may be 
evaluated in any number of dimensions. Consequently, the value of the RBF is not 
dependent on direction. 

There are numerous radial basis functions in common use for a range of 
applications. Some of the most common are Gaussian, multi-quadric, and exponential. 
These functions are generally chosen to exhibit some smooth, continuously decreasing, 
variation in value as the distance from the center increases. For interpolation 
problems, it is also desirable that the functions do not have a compact support, i.e. 
the radial basis function should have a non-zero value across the entire interpolation 
region. When RBF have infinite support, all RBFs extend through all data sites and 
interpolants. The weights of the RBF can then be determined by solving a linear 
system, where the vector of weights is solved for from a coefficient matrix of RBF 
values for each combination of center and data site location and the vector of known 
function values, /: 


<F(r(a;i,ci)) <F(r(a;i, C 2 )) ... <F(r(a;i, cat)) 
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$(r(xAr,ci)) $(r(a;7v,C2)) ... <h(r(a;Ar, ctv)) 




f{xN) 


where, for example, the multi-quadric RBF (shown in Figure [3?f| ) has the form: 

<F(r(a;, c)) = a/ 1 -|- (er(a;, c))^ (3.6) 

This formulation produces a function approximation which is exact at all known 
data points. If the RBF are chosen with compact support, some data sites may not 
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Multi-Quadric Radial Basis Function 



Figure 3.7: The multi-quadric Radial Basis Function (RBF). 


lie within some RBF and the resulting problem may be over or under constrained. 
Additionally, approximation with compactly supported RBF limits the domain of 
the approximating function since the RBF have non-zero values over a hnite region. 

The centers of the radial basis functions may be chosen arbitrarily; however, 
for interpolation problems good results are generally achieved by collocating the 
RBF centers with the data site locations. This collocation ensures that each data 
site is strongly represented by at least one RBF. The formulation of the RBF-based 
approximation does not depend on the distribution or uniformity of data sites. An 
example of the multi- quadric RBF used to approximate the sine function using hve 
bases is shown in Figure |3^ 
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Figure 3.8: The approximation of the sine fnnction nsing Radial Basis 
Functions (RBF). The black line is the function sin27ra;, the dashed lines 
are the hve RBF bases, and the red line is the RBF approximation using 
the superposition of the bases. 
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3.3.3 Neural Network Estimation 


An Artificial Neural Network (ANN) [21] can be seen as a more general version 
of the RBF function approximation discussed in the previous section. In the context 
of an ANN, the basis functions are called neurons. Each neuron accepts a set of 


parameters as an input, and produces one output value, shown in Figure [3i9] — much 
like a single RBF. In the RBF function approximation, the weighted sum of these 
outputs is taken as the overall output of the RBF approximation; this is equivalent to 
an ANN with a single layer of neurons, as seen in Figure |3A0 However, in a typical 
ANN there are multiple layers of neurons. The hrst layer accepts the inputs and is 
called the input layer. The last layer of yields the output values and is therefore called 
the output layer. The intermediate layers of neurons are called hidden layers — these 
layers accept the output of another layer as input, and likewise pass their output 
to another layer as inputs. Weights are determined for every input-output pair in 
the network, providing the ANN great flexibility in modeling nonlinear phenomena. 
Generally, the neurons can be described by any multi-input single-output function, 
although RBF are commonly used because of their utility in function approximation. 


A diagram of a simple ANN with two hidden layers is shown in Figure 3.11 


The process of determining the best ANN weights in order to £t the behavior 
of the ANN to known data is called training. There are a number of schemes in 
use to train ANN, most of which apply optimization techniques like those discussed 


in Section 3^ in order to minimize the error between the ANN output and known 
training data for a given set of inputs. This error is dehned as the mean sum of the 
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Output 


Figure 3.9: A single neuron. 


squares of the error at each of the data points in the training set, i.e.: 

1 ^ 

F '^{fANwiXi) - ftrain(Xi)Y (3.7) 

i=l 

However, due to the flexibility of ANN they are susceptible to “over-htting” , 
i.e. the ANN accurately reproduces the training data, but provides inaccurate results 
for inputs outside the set used in the initial training. The most common technique 
to improve the generalization of ANN to unknown data is called early stopping. IS] 
In this process, some randomly selected portion of the known data is withheld from 
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Figure 3.10: A single-layer Radial Basis Function (RBF) network. 


the training set, and is called the validation set. On each iteration of the training 
process, the output of the ANN is compared against the validation set — initially, the 
error on the validation set decreases along with the error on the training set. At 
some point in the process, the error on the validation set begins to increase while 
the error on the training set continues to decrease. This indicates that over-htting is 
beginning to occur, and the training process is then terminated. Early stopping is an 
effective technique for improving the generalization of ANN, but is not well suited 
to limited data sets, since some data must be set aside for validation. A different 
approach which results in slower training, but better generalization for smaller data 
sets, is called regularization. Regularization changes the metric used to assess the 
error of the ANN during the training process, including a term which limits the 
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Figure 3.11: An Artificial Neural Network (ANN) network with two 
hidden layers. Fully connected network relations are simplified with block 
arrows. 


magnitude of the network weights: 

N n 

^ ^ “ ftrain{Xi)f + (1 ~ (3-8) 

i=l ^ j=l 

f3 is called the performance ratio. Determining the optimal value of for 
regularization a priori is difficult, but techniques exist in order to automatically 
adapt the performance ratio in order to ensure effective regularization. One of 
the most effective and commonly employed automatic regularization techniques is 
Bayesian regularization, |53] which uses probability theory and statistical arguments 
to estimate how much regularization is warranted based on the sampling of data in 
the training set. In this dissertation, the Bayesian regularized Levenberg-Marquardt 
training algorithm provided in the MATLAB® Neural Network toolbox[5l] is utilized 
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to fit the ANN which relate the governing parameter valnes to the dependent 
parameter values discovered from the measured data using parameter identihcation 
technique described earlier in this chapter. The combination of the trained ANN and 
the analytical model which accepts the dependent parameters results in a complete 
and generalizable model of the noise radiated by a particular rotor conhguration 
under a wide range of non-dimensionally defined operating conditions. 

3.4 Model Usage 

Once suitable models have been identihed describing the external noise radiation 
of the helicopter, they can be applied to estimate noise at any specihed flight condition. 
First, the specified flight condition must be defined in terms of the non-dimensional 
governing parameters of each noise source, e.g. main rotor harmonic noise, tail rotor 
harmonic noise. These parameters can be readily determined from the dimensional 
parameters dehning the flight trajectory and ambient conditions to be used in 
the ground noise contour simulation. The non-dimensional governing parameter 
values are then used to select dependent parameter values for the assumed models 
using the parameter estimation scheme describing the mapping between governing 
and dependent parameters. If the multi-objective parameter identihcation scheme 
provides a Pareto frontier of possible solutions, the set of dependent parameter values 
can be chosen at runtime based on the type of helicopter operation performed. The 
noise radiated to ground based observers is then computed for each noise source 
using the set of governing and dependent parameters for each noise source. The 
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output of each noise source model is superimposed to yield a complete description of 
noise radiation at the desired operating state. 
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Chapter 4 


Measurement and Data Processing 

4.1 Source Noise Measurements 

The simplest way to gather noise measurements from a wide range of directions 
for a helicopter in a steady flight condition is to utilize an array of ground-based 
microphones. For straight-line flight conditions, a linear array of microphones has 
proven effective, as is typically employed to collect data for use with RNM. As the 
helicopter passes over the linear array, shown in Figure |4T| the microphones measure 
noise radiated from the front of the helicopter while the vehicle is ahead of the array 
and aft of the helicopter after it passes the array. These measurements are made 
across a range of lateral directivity angles by positioning a number of microphones 
along the ground in a line orthogonal to the flight path of the vehicle. The precise 
geometry of the measurement can be determined as a function of time given position 
tracking data from the helicopter time-synchronized with the acoustic measurement 
as well as the locations of the stationary microphones. 

4.2 Virtual Inflight Observers 

Measurements of noise on the ground are not easily comparable. The distance 
between the source and the observer varies with directivity angle, as well as with 
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Figure 4.1: The linear microphone array used to construct acoustic 
hemispheres for Rotorcraft Noise Model (RNM). 


changes in the trajectory of the aircraft. Differences in the propagation distance result 
in differences in observed noise levels due to spherical spreading and atmospheric 
absorption. In addition, the relative motion between the aircraft and the observers 
resnlts in an apparent freqnency shift in the observed aconstic signals. In order 
to have directly comparable sets of measnrements, the measnrements shonld be 
transformed to a common frame of reference. 

The acoustic measurements can be transformed to a common frame of reference 
by introdncing the concept of virtnal inflight observers. jSS] The virtnal inflight 
observer makes the assumption that the helicopter can be represented as a compact 
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Figure 4.2: Doppler effect on acoustic pressure time histories measured 
by stationary observer. 


noise source, with all noise from one rotor radiating from a single point on the rotor 
disk. The virtual inflight observer then travels through the medium with the source 
maintaining hxed distance from the source at all times. This is similar to a true 
inflight observer who is flying alongside the helicopter and analogons to measurements 
of a hxed source from hxed observer locations in the presence of a moving medium, 
e.g. acoustic measurements made in a wind tunnel. However, because the virtual 
inhight observer is composed from the transformation of ground-based acoustic 
signals, and the helicopter noise source is highly directional, the virtual inhight 
observers must orient themselves with respect to the helicopter noise source so as to 
match the corresponding directions between the noise source and the ground-based 
observers. 

The transformation of the acoustic pressure time-history signals from from 
ground-based microphone measurements to those measured by the virtual inhight 
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Figure 4.3: Diagram of the virtual inflight observer concept. 


observers is accomplished through a time-domain process of de-Dopplerization. Unlike 
traditional de-Dopplerization approaches, which are formulated in the frequency 
domain, inn] the time-domain approach fully models how the acoustic signal deposited 
in the medium is perceived due to the relative motion between the source and observer. 
This process does not attempt to compensate for the effect of Doppler amplification, 
which is due to the motion of acoustic sources with respect to the medium itself. 
Without a complete understanding of the type, distribution, and motion of the 
individual sources through the medium, it is impossible to accurately model the 
Doppler amplihcation process. However, this process is fully accounted for in the 
physics-based analytical modeling structure described in Chapter during the 
evaluation of the Ffowcs Williams - Hawkings equation, as an estimate of the 
composition of sources can be made in this case using the parameter identihcation 
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techniques of Chapter 


The time-domain de-Dopplerization process is a simple accounting of the 
retarded time equation, which describes the time delay between the emission of noise 
into the medium and reception from the medium by the observer. 


t — T — 


Oq 


0 


(4.1) 


where t is the time of reception, r the time of emission, and \x — Xs\ the distance 
between the source and observer through the medium. The theoretical development 
of this transformation can be found in Appendix [Aj 



Figure 4.4: The process of time-domain de-Dopplerization. 


In time-domain de-Dopplerization, this equation is evaluated continuously over 
the entire length of the recorded signal. This requires an accurate knowledge of the 
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position of the helicopter as a function of time. The process of de-Dopplerizating 
the pressure time-history data is as follows: 

• The helicopter pressure time-history and position tracking data are discretized 
in time at hxed sampling rates; in general, the position tracking data is 
discretized at a much lower sampling rate. 

• In order to provide tracking data at arbitrary time steps a continuous piecewise- 
dehned spline is £t to the tracking data. 

• For each acoustic sample, the correct time of emission, time of observation 
at the spherical surface, and position of the helicopter along the track are 
calculated. For a real trajectory, this solution is iterative but convergence 
is quick. The iterative solution may be initialized by assuming straight-line 
motion of the helicopter between known position tracking points, for which 
the retarded time is known exactly. 

• Using the propagation distance calculated at each time-history sample, the 
change in amplitude due to spherical spreading may be calculated and applied 
for each individual pressure time-history sample. The de-Dopplerized signal 
observed by the “virtual inflight microphone” may now be obtained. 

• The time steps between samples of the transformed signal are non-uniform. In 
order to calculate frequency spectra using standard techniques the signal must 
be resampled to a hxed time-step. This process may introduce aliasing at a 
lower frequency than for the untransformed signal. An increase of the initial 
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sampling rate by 30% has been recommended. However, the Non-Uniform Time 
Discrete Fourier Transform [37] can be used to calculate uniform frequency 
spectra without affecting the Nyquist frequency, as the average sampling 
frequency of the non-uniformly spaced signal does not change significantly. 

The result of this process is pressure time-history signals as observed by virtual 
inflight microphones which slowly trace paths along the surface of a hemisphere set 
a hxed distance away from a point on the rotor disk, corresponding to the directions 
measured during the data gathering run. Assuming that the helicopter maintained 
a steady flight condition throughout the run, this forms a complete representation 
of the noise radiated by the helicopter in terms of direction, frequency, amplitude, 
and phase during that flight condition. However, in order to apply the parameter 
identihcation technique of Chapter to these data, the contributions of the various 
rotor harmonic noise sources must be separated from the transformed pressure 
time-history signal. 
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4.3 Rotor Harmonic Noise Separation 


A common method for separating rotor harmonic noise from measurements 
made in a wind tunnel is periodic ensemble averaging. [58] This method has also 
been applied to inflight measurements of the noise radiation of entire rotor craft. [5^ 
The acoustic pressure time-history signal is divided into individual rotor revolutions. 
Since the rotor rotational rate is not necessarily constant, data from a once per 
revolution sensor mounted to the rotor shaft is generally utilized to identify the 
actual rotor periods. Under the assumption that the acoustic state of the rotor is 
generally stationary from one revolution to the next, the pressure time history signals 
from a number of rotor revolutions are scaled to a normalized time and ensemble 
averaged. The averaged signal retains the component which is harmonic to the rotor, 
but suppresses the part of the signal which is not harmonic to the rotor. The number 
of rotor revolutions included in the ensemble average determines the quality of the 
suppression of non-harmonic noise. For wind tunnel measurements, typically on the 
order of fOO rotor revolutions are used to generate the separated signal. |5H| The 
process can be applied to both the main and tail rotor to identify the contribution 
of each to the overall external noise radiation. Subtracting the two signals from the 
original measured signal produces a residual pressure time-history that includes the 
contribution of all non-rotor harmonic noise sources. 

Having introduced the technique for transforming ground-based acoustic pres- 
sure time-history signals to those measured by virtual inflight observers, it is possible 
to perform periodic ensemble averaging on measured ground-noise data in order 
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Figure 4.5: Periodic ensemble averaging of rotor noise in order to suppress 
the contribution of non-rotor harmonic sources. 


to separate the contributions of main, tail, and non-rotor harmonic noise sources 
to the overall noise radiation, as would be performed in a wind tunnel. [33] The 
individual rotor revolutions must hrst be identihed in the transformed signal, so as 
to compensate for variations in the rotor rotational speed. 

For routine measurement of helicopter noise, a once per revolution sensor is 
not always available to provide direct measurements of the rotor periods. In this 
case, a virtual once per revolution sensor can be constructed from the transformed 
pressure time-history acoustic signals. This is accomplished by calculating a wavelet 
power spectrogram of the transformed signal. The continuous wavelet transform is 
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expressed using the following convolution: 



t-t' 


a 


dt' 


(4.2) 


where a is the scale of the wavelet (inversely related to frequency), t is the time 
of the transform, and is the complex conjugate of the wavelet basis function. In 
this dissertation, the Ricker wavelet basis function is used, as it is known to provide 
excellent temporal accuracy. [60] This wavelet basis is expressed as: 


where a is a shape parameter, typically chosen about unity. The wavelet power 
spectrogram differs from traditional short time Fourier transform techniques in that 
an arbitrary frequency resolution can be employed to analyze the signal. For this 
application, the spectrogram is calculated at a number of wavelet scales associated 
with the frequencies between 90% and 110% — by picking out the maximum power 
over this frequency interval in time, the effect of individual blade passages can be 
identihed, yielding the corresponding rotor period intervals. The entire rotor period 
identihcation process is shown in Figure [4.6 on the next page 

Rotor periods are collected for averaging based on angle between the prop- 
agation path from the assumed source and observer and the rotor tip-path-plane. 
Since the virtual observers are composed from ground-based measurements, they 
traverse through a range of elevation angles over time. At positions near the horizon, 
the rate of angular motion is small, since the ground-based microphones were far 
from the helicopter during the measurement. At positions underneath the helicopter. 



(4.3) 


the angular rate of change of the virtual observers is more rapid, as the signal was 


De-Dopplerized Pressure Time History Signal 



0.5 1.0 1.5 

Wavelet Power Spectrogram 



Rotor Period Detection Signal 



Figure 4.6: Calculation of the wavelet power spectrogram to identify 
rotor blade passages from the de-Dopplerized acoustic signal. 
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collected as the helicopter was passing directly over the ground-based microphone 
array. In order to correctly reflect noise radiated from a single direction, the periodic 
ensemble averaging process must be limited to a small angle swept by the virtual 
observers. Representative main and tail rotor pressure time history signals are 
created using periodic ensemble averaging at fixed angular increments along the 
paths traced by the virtual observers. Consequently, more rotor revolutions may be 
included in the average for virtual observers oriented towards the horizon than while 
underneath the helicopter. (Figure |47|) 


Noise Source at 
Time of Emission 



Figure 4.7: Diagram illustrating the the inclusion of rotor periods in the 
ensemble averaging process as a function of observation angle. 


So far, the transformation process has assumed that all noise produced by the 
helicopter may be represented by a compact source located, for instance, at the main 
rotor hub. In reality, the noise is generated by a distribution of sources across the 


rotor blades which vary azimuthally in strength. The assnmption of a hnb-centered 
compact sonrce introdnces an error in the phase of the aconstic pressnre originating 
from sonrces away from the hnb dnring the evalnation of the retarded time equation. 
For lower frequency sources, such as steady loading and thickness noise, the error in 
phase is relatively small. However, for impulsive noise sources like those associated 
with BVI, even a small change in phase can disrupt the averaging process. The error 
in phase is dependent on the geometry between the observer, the assumed source 
location, and the real source. 



Figure 4.8: The phase error due to the compact source assumption 
changes with respect to observation angle. 


Figure |T8 shows the how the compact source assumption introduces a phase 


error during the transformation process due to the longitudinal distance between 
the assumed source location and any individual real source on the rotor. When the 
observer lies along the line containing the assumed source location and the real source, 
the error in retarded time calculation will be the distance between the two divided 
by the speed of sound. When the observer is directly underneath the two sources, 
the error in the retarded time calculation will be zero, as the path lengths will be 
equal. If the observer is much farther from the sources than the distance between 

the real and assumed source, the error in retarded time may be approximated as: 

At =— cos 6 (4.4) 

Oo 


where d is the longitudinal distance between the assumed and actual source 
location, and 6 the angle between the observer and the rotor plane. Plotted against 


elevation angle, this results in the phase error seen in Figure 4.9 on the next page 


More significantly for the averaging process, the rate of change of the error in 
the retarded time calculation is highest out of the plane between the assumed source 


and the real source, shown in Figure 4.10 on page 92 


For the case of a real helicopter, all sources are on or near the rotor blades. 
Several different BVIs may occur at once, originating from different locations on the 
rotor disk. At the elevation angles underneath the rotor plane, the phase of the BVI 
pulse will shift quickly between each rotor period, eliminating the BVI pulse from 
the averaged signal. The error in the assumed source location is not as important 
in the plane of the rotor because the rate of change of phase due to error in the 
retarded time calculation is small. Therefore, the BVIs which radiate the most noise 
below the helicopter are the most sensitive to error in source location. The shifting 
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Figure 4.9: Error in time with elevation angle due to compact source assumption. 


phase of such a BVI from one revolution to the next is shown in Figure |4.11 


The hemisphere could be constructed about a different assumed source location 
on the rotor, more representative of BVI. However, different types of BVI occur 
during different flight conditions, occurring at different areas on the rotor disk. For 
example, a parallel BVI may occur on a 2-bladed rotor at about 45° azimuth, and 
occurs across much of the blade simultaneously. Oblique BVIs occur across a range 
of azimuth, and sweep across the blade over time. Since it is difficult to predict 
which BVI will be dominant for a particular flight condition a priori, it is useful 


91 



Figure 4.10: Rate of change of time error with respect to elevation angle. 


to correct for the error in the phase of BVI pulses individually at each directivity 
angle on the snrface of the aconstic hemisphere. The most prominent BVI pnlses 
in the signals are synchronized between revolutions by adjusting the phase of each 
time segment nsed in the periodic average so as to align the peak response of the 
cross correlations between the signal for each rotor period and the signal of the 
central rotor period in time. The result is a representative rotor harmonic pressure 
time-history signal containing both the impulsive and lower harmonic noise radiated 


in the corresponding direction, shown in Fignre 4.12 on page 94 This process is 
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Figure 4.11: Changing time error causes impulsive BVI noise to shift in 
phase over successive rotor time periods. 


repeated for each direction traversed by the virtual observers in order to construct 
a complete representation of rotor harmonic noise radiated by the helicopter. By 
repeating the process for both the main and tail rotor, all rotor harmonic noise 
sources can be characterized. Non-rotor harmonic noise is captured by subtracting 
the rotor harmonic noise signals for each rotor from the imaveraged virtual observer 
pressure time-history signals in order to yield a residual signal. This residual signal 
can be evaluated using standard frequency- domain techniques for generating acoustic 
hemispheres, since phase does not carry useful information for these noise sources. 
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Phase Corrected Time History Window 



Figure 4.12: A cross-correlation can be used to adjust the phase of each 
rotor period so as to preserve the BVI impulse in the averaged signal. 
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4.4 Projection of Acoustic Hemispheres 


Three-dimensional images of hemispheres are unable to show the noise levels 
across the entire surface of the sphere without employing multiple views. In order to 
display the hemisphere in its entirety on a single plot, it is necessary to project the 
hemisphere onto a 2D planar surface. 

No 2D projection of a spherical surface can be free of distortion, as Euler proved 
in 1777.(61] Such a projection can either be conformal or equal-area, but never both. 
A conformal projection preserves the shape of features, since at each point the angles 
between the spherical surface coordinates are locally preserved in the flat coordinate 
system. However, the relative size of features will be changed globally during the 
transformation. An equal-area projection ensures that each region of the map will 
retain the same area in the flat coordinate system as it had on the spherical surface; 
however, the orientation of points with respect to each other within any region will 
be changed, distorting the shape of each feature. In order to effectively display all 
of the information captured by a hemisphere, some sort of projection to a flat 2D 
figure must be employed. 

In this paper, hemispheres are displayed mapped to a 2D image using a Lambert 
Conformal Conic projection, developed and used commonly by cartographers to 
accurately display maps of polar regions. The Lambert projection is perfectly 
conformal, but is also fairly close to being equal-area when no more than one half 
of a sphere is transformed at once. This is achieved by splitting the back end 
of the hemisphere along the 0 degree azimuth angle, and unrolling the resulting 
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cone to a 2D projection. (Fignre 4.13) The top of the projection represents the 
angle directly in front of the helicopter, the split seam the angle directly behind, 
the right side the advancing side of the helicopter, and the left the retreating side. 
The center of the projection, at 90 degrees elevation, represents the bottom of the 
hemisphere. The radial lines thus represent the azimuth angles, and the concentric 
lines the elevation angles. This projection allows for an intuitive understanding of 
the directivity patterns of externally radiated noise. 


96 


Front ^=180*, ^=0* 




0 * 


Figure 4.13: The Lambert projection of an acoustic hemisphere. 
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4.5 Measurement Quality 


One fundamental assumption of all processes for characterizing helicopter noise 
from ground-based measurements is that the helicopter maintains a steady acoustic 
state over the measurement window, i.e. the noise radiation of the helicopter is 
identical from one rotor revolution to the next. In practice, this is never the case — the 
flight state of the helicopter is continuously changing due to atmospheric unsteadiness 
and pilot control corrections. It is important to ensure that the consequent variability 
in the acoustic state of the helicopter is within acceptable limits throughout the 
measurement run. Without simultaneous measurement of the helicopter at fixed 
observation angles, it is difficult to directly quantify the quality of the measured 
data. However, some indirect observations may be made to place bounds on the 
acoustic variability. 

One simple approach is to repeat each measurement multiple times, and assess 
the variability of the measurement statistically. Provided that similar noise levels 
are measured over all directions for several repeated runs, the standard deviation 
in noise levels can be calculated and a level of confidence can be established in the 
quality of the measurements, for instance by calculating the standard deviation of 
the SPL at each point on the resulting acoustic hemisphere. The degree of confidence 
increases as more measurements are included in the analysis. However, this requires 
that each flight condition to be characterized must be measured many times, greatly 
increasing the time and expense of collecting helicopter noise data. 

Another statistical method of characterizing helicopter noise is to take a moving 



standard deviation of the noise levels seen by the virtual observers as they traverse 
the acoustic hemisphere surrounding the helicopter. With the knowledge that noise 
levels do not change sharply with directivity, it is assumed that sudden changes in 
noise levels found in the transformed signals are due to changes in the flight state of 
the helicopter, and can be captured by taking calculating the standard deviation 
over a short time window. This method provides some estimate of the variability of 
the acoustic state of the helicopter from a single measurement, and also indicates 
the quality of the measured data as a function of directivity angle. The 2N point 
moving standard deviation at point k is formulated as: 

^ k+N 

= \ ^ - XnY (4.5) 

\ n=k-N 

where Xn is the central moving average at point n. 

Both of these techniques will be applied to assess the quality of the measured 
data in the following chapter. 
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Chapter 5 


Applications 

The data analysis and modeling methods develop in Chapters and can be 
applied not only to ground-based measurements of flying vehicles, but also to wind 
tunnel test data for isolated rotors. In this chapter, models are built using two sets of 
data for two-bladed rotors — ^the first is from the test of an isolated 1/7^^ scale model 
of the AH-1 Operational Loads Survey (OLS) rotor in the German-Dutch Windtunnel 
(DNW).[63] The second is from the flight testing of a Bell 206B3 helicopter over 
Moffett Field, California. |M] 


5.1 Wind Tunnel Test 


The parameter identihcation method is hrst applied to acoustic measurements 
of several flight conditions produced from the test of the 2-bladed 1/7^^ scale AH-1 
Operational Loads Survey (OLS) rotor in the German-Dutch Wind Tunnel (DNW). 
This dataset is notable for high quality acoustic measurements over a wide range 
of operational conditions carefully controlled by the governing non-dimensional 
parameters of rotor harmonic noise sources. Nine microphones located within the 
wind tunnel jet, labeled 2 through 10 in Figure |5.1 on page 105 , were included 
in the parameter identihcation process. Although these microphone locations are 
predominantly ahead of and towards the advancing side of the rotor, there is still a 
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R 

3.14 ft 

c 

0.341 ft 

t/c 

9.71 % 


Table 5.1: Geometric characteristics of the model scale OLS rotor. 


sufficient variation in directivity angles to allow a single “best ht” to be identihed 
in all cases. Using these data, the ability of the parameter identihcation techniqne 
to produce models correctly capturing the parametric variations of rotor noise is 
evaluated. The model scale OLS rotor has the geometric characteristics listed in 
Table 15.11 

The data collected during this test is very consistent, especially for the mi- 


crophones inside of the wind tunnel jet used in this research. Figure 5.2 compares 


typical 100 revolution averaged and instantaneous pressure time-history measure- 
ments recorded from microphone channel ^4 dnring the experiment. There is little 
difference between the instantaneous and averaged measurements, especially in terms 
of the impnlsive BVI noise — this demonstrates the steadiness of the rotor acoustic 
state during the test. This steadiness is one of the major benehts of a carefully 
controlled wind tunnel test. The baseline operating condition is dehned by the 


non-dimensional governing parameter valnes in Table 5.2 


Figures 5.3, 5.4, and 5.5 show the measured and simulated acoustic pressure 


time-histories for the baseline case and htted model at in-plane, 30° and 45° ont-of- 
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A 

0.0135 

Ct 

0.0054 


0.164 

Mht 

0.620 


Table 5.2: Non-dimensional governing parameter values for wind tunnel test cases. 


plane observer locations, respectively. In general, the pressure time-histories for most 
observer locations are fitted by the model to within about 1 dB BVISPL; although 
some of the in-plane observers are ht to within 1.5 dB BVISPL. These in-plane 
observers capture the lower harmonic thickness noise well. In addition, the overall 
magnitude of the in-plane BVI noise is fairly well captured by the htted model, but 
the phasing of the impulse is not in complete agreement with the measured data. 
The more out-of-plane observers accurately capture form of the BVI impulses in the 
pressure time history signal. Lower harmonic noise is well captured by the model as 
well, with OASPL also ht to within 1 dB (not shown) for all observers. Overall, the 
discrepancy between the measured and modeled pressure time-history signals is well 
within the natural variability of noise from one rotor revolution to the next observed 
in actual hight. 

As an accurate analytical model of the external noise radiation of the rotor 
for a single hight condition has been developed, the generalization of this model to 
other operating conditions is now demonstrated. Using the modeling parameters 
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identified for the baseline condition of the model rotor estimates are made for the 


noise radiated at other operating conditions, as dehned by variations of the non- 
dimensional governing parameters. These noise estimates are then compared to 
measured acoustic data for the OLS rotor undergoing the same parametric variations 


of operating condition. Figure 5.6 shows the peak-to-peak acoustic pressure for 


various thrust coefficients as seen at the 45° out-of-plane observer location shown 
for the baseline case in Figure |53j Both the measured and simulated data show an 
approximately linear trend in peak levels with variation in thrust coefficient; however, 
the slope of the simulated data is signihcantly less than that of the measured data. 
On closer inspection, the linear ht of the simulated model extrapolates to zero peak 
acoustic pressure at zero thrust, unlike the linear ht of the measured data. This 
indicates that some physics of the BVI noise generation process are not captured by 
the analytical model. 

Similarly, the peak acoustic pressure trend for the measured and modeled data 
due to variations in advancing tip Mach number is shown in Figure |T7 The model 
has a similar trend to the measured data, but over-predicts noise at lower Mach 
numbers and under-predicts at higher Mach numbers than the baseline case. This 
is because the chosen analytical model does not directly incorporate the effects of 
compressibility on BVI. 

In order to model physics not explicitly captured in the analytical model, the 
dependent modeling parameters must be varied as a function of the non-dimensional 
governing parameters. This can be accomplished by identifying dependent parameters 
for several different operating conditions and then using this information to estimate 
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the dependent parameter valnes for other nnmeasnred operating conditions — this 


process is described in Section 3.3 This resnlts in a large improvement in the 


accuracy of the noise models when generalized to unmeasured operating conditions. 


For example, Figure 5.8 shows the results of using simple linear interpolation to 


estimate dependent parameter values from measured data for two different different 
thrust coefficient values. The discrepancy between measured and modeled peak 
acoustic pressure values is greatly reduced over those observed when using only the 


dependent parameters developed from the baseline condition, as shown in Figure 5.6 


Adding just one more additional set of identihed parameters further improves the 
ht, allowing the model to match the observed trend almost exactly, as shown in 


Figure |5.9 

Likewise, incorporating additional measured data can also improve general- 


ization of the model as a function of advancing tip Mach number. Figures 5.10 


and 5.11 show the peak BVI acoustic pressure variation as a function of advancing tip 


Mach number for models incorporating data from two and three different operating 
conditions, respectively. Model agreement with measured data improves dramatically 
by incorporating just one additional set of measured data; further rehnement is 
achieved by adding data for additional Mach number variations. 
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Figure 5.1: Side-view of the DNW windtunnel microphone measurement 
configuration for the 1/7^^ scale OLS rotor. 
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Acoustic Pressure, Pa 


DNW 1/7 Scale OLS Pressure Time History Comaparison 



Figure 5.2: Comparison of typical averaged and instantaneous pressure 
time-history measurements for the DNW 1/7^^ scale OLS experiment. 
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Acoustic Pressure, Pa 


DNW 1/7 Scale OLS Microphone Channel #2 Fit, ABVISPL = -1 .348 



Figure 5.3: Measured and fitted-model acoustic pressure time-histories 
for in-plane observer location of rotor in baseline operating condition. 
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Acoustic Pressure, Pa 


DNW 1/7 Scale OLS Microphone Channel #3 Fit, ABVISPL = -1 .006 



Figure 5.4: Measured and fitted-model acoustic pressure time-histories for 
30° out-of-plane observer location of rotor in baseline operating condition. 
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Acoustic Pressure, Pa 


DNW Ml Scale OLS Microphone Channel #4 Fit, ABVISPL = -0.2863 



Figure 5.5: Measured and fitted-model acoustic pressure time-histories for 
45° out-of-plane observer location of rotor in baseline operating condition. 
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Peak Acoustic Pressure, Pa 



Thrust Coefficient 


Figure 5.6: Measured and simulated variation in peak BVI acoustic 
pressure with thrust coefficient for the 45° out-of-plane observer. Model 
uses parameters developed for the baseline case. 
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Peak Acoustic Pressure, Pa 



Advancing Tip Mach Number 


Figure 5.7: Measured and simulated variation in peak BVI acoustic pres- 
sure with advancing tip Mach number for the 45° out-of-plane observer. 
Model uses parameters developed for the baseline case. 
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Peak Acoustic Pressure, Pa 



Thrust Coefficient 


Figure 5.8: Measured and simulated variation in peak BVI acoustic 
pressure with thrust coefficient for the 45° out-of-plane observer. Model 
uses parameters developed for the baseline case and Ct = 0.0080. 
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Peak Acoustic Pressure, Pa 



Figure 5.9; Measured and simulated variation in peak BVI acoustic 
pressure with thrust coefficient for the 45° out-of-plane observer. Model 
uses parameters developed for the baseline case, Ct = 0.0047 and Ct = 
0.0080. 
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Peak Acoustic Pressure, Pa 



Advancing Tip Mach Number 


Figure 5.10: Measured and simulated variation in peak BVI acoustic pres- 
sure with advancing tip Mach number for the 45° out-of-plane observer. 
Model uses parameters developed for the baseline case and Mat = 0.64. 


114 


Peak Acoustic Pressure, Pa 



Advancing Tip Mach Number 


Figure 5.11: Measured and simulated variation in peak BVI acoustic 
pressure with advancing tip Mach number for the 45° out-of-plane ob- 
server. Model uses parameters developed for the baseline case, Mat = 
0.64 and Mat = 0.84. 
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In addition to producing noise estimates at known observer locations, the 


models developed by this process can be used to estimate noise levels in directions 


which were not measured. Figure 5.12 shows the Lambert projection of the BVISPL 


noise level contours produced using the model for the baseline case on the surface of 
a hemisphere offset 30 feet from the rotor hub. In the flight condition, BVI noise is 
radiated ahead of and below the rotor disk. The Lambert projection is described in 


more detail in Section 4.4 As advance ratio increases, (Figures 5.13 through 5.15) 
the noise levels increase and are directed more in-plane and towards the advancing 


side, as has been observed in previous theoretical |65j and experimental studies of 


BVI.[5B] 
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Figure 5.12: BVISPL noise level contours for the baseline case (/r = 0.164). 
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Figure 5.13: BVISPL noise level contours for /r = 0.194. 


118 




119 



Figure 5.15: 


BVISPL noise level contours for /r 


0.270. 
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5.2 Flight Test 


In addition to wind tunnel data, the modeling approach is applied to measured 
ground noise data from the flight test of a Bell 206B3 helicopter at Moffett Field, CA, 
In 2006. Ground noise data were collected from a linear array of eight ground board 
microphones for a variety of steady and maneuvering flight conditions. In addition, the 
helicopter was Instrumented with a unique Inflight measurement array. The modeling 
approach developed In this dissertation Is applied to this dataset to demonstrate how 
analytical models that relate noise to non-dlmenslonal governing parameters may be 
built for full-scale helicopters. Such models are suitable for ground noise contour 
generation under a wide range of flight conditions. Ground noise contours for a 
diverse range of atmospheric and flight conditions are often required for effective 
mission planning and community noise abatement, yet measured ground noise data 
are generally available over only a limited set of operating conditions due to the 
significant cost of acquiring such data. Non-dimensional analytical models provide 
a means to generalize a limited set of measured data so that it is applicable over 
a wider range of conditions. In addition, the relatively simple analytical models 
resulting from this technique can be run in real time, which makes them suitable for 
Inflight noise estimation. The Bell 206B3 main rotor has the geometric characteristics 
listed in Table [T3l 

As described In the previous section, the ground-based acoustic measurements 
are transformed to a reference frame analogous to wind tunnel measurements. The 
wavelet transform is then used to identify individual rotor revolutions from the 
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R 

16.7 ft 

c 

1.10 ft 

tjc 

11.3% 


Table 5.3: Geometric characteristics of the Bell 206B3 main rotor. 


acoustic signal. The main rotor harmonic pressure time-history signals are extracted 
using a periodic averaging process for a number of observer directions. More rotor 
periods may be included in the averaging process for in-plane observers than out-of- 
plane observers; therefore, the more in-plane data is a better representation of the 
main rotor harmonic time-history signal. In addition, the dominant noise sources 
of the rotor do not strongly radiate directly underneath the helicopter. For these 
reasons, the distribution of “virtual inflight observers” included in the analysis is 
biased away from the underside of the helicopter and towards the dominant radiation 
directions. The measured signals transformed to these observers are the inputs to 
the parameter identihcation process. 

The modeling process is applied to data for several steady descending flight 
conditions at 60 kts airspeed with approximately the same gross weight and rotor 
RPM for each case — in effect, a parametric variation of non-dimensional inflow which 
greatly affects BVI noise levels. The non-dimensional governing parameter values for 
this case are shown in Table 15.41 

Due to the increased variability in measured flight test data, the model does 
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A 

various 

Ct 

0.0033 

/i 

0.15 

Met 

0.626 


Table 5.4: Non-dimensional governing parameter values for flight test cases. 


not match the measured data as closely as for the wind tunnel case. However, the fit 
is within 1.5 dB BVISPL for all observer locations — well within the typical variations 
between repeated measurements of the same operating condition in flight test data. 


Figure 5.16 shows the unaveraged BVISPL levels for the centerline “virtual observer’" 


at 5° increments for three different runs of the nominal 60 kts, -6° flight path angle 
flight condition. The condition was held fairly close to the desired test point in all 
three cases. In Run A, the average airspeed recorded was 61 kts. In Run B and Run 
C, the average airspeed was approximately 56 kts. The average flight path angle in all 


cases was within 0.2° of nominal. Figure 5.17 shows a plot of the standard deviation 


between all three runs at each observation angle. Significant deviations in measured 
BVISPL levels can be observed for all three cases in the range of 1. 0-2.0 dB BVISPL, 
especially near the front of the helicopter, where the BVI noise levels are the highest. 
This is about the same level of repeatability observed for approach conditions in other 


rotorcraft flight test measurement programs. (6^ Figure 5.18 shows the hve-point 
moving standard deviation of the same “virtual observer” measurement — it can be 
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seen from this calculation that the largest variability in noise radiation over time time 
seems to occur underneath the helicopter, where retreating side BVI might contribute 
to the BVISPL noise levels. Of these three repeated cases, Run A was selected for 
modeling, due to low variability observed in the moving standard deviation and a 
higher “steadiness” of the condition throughout the run, as indicated by the author 
and the flight director from the helicopter during the run. 



Observation Angle Below Horizon, degrees 


Figure 5.16: Unaveraged BVISPL at the centerline “virtual observer” 
during the 60kts, -6° flight condition. 


Hemispheres developed by averaging and interpolating the measured main 
rotor harmonic noise in the BVISPL frequency range are shown for the -6° and 
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Figure 5.17: Standard deviation of BVISPL between Runs A, B, and C 
at the centerline “virtual observer” during the GOkts, -6° flight condition. 


-7.5° flight path angle cases in Figure and Figure [O0| respectively. Due to the 
high drag of the helicopter in the experimental conhguration, BVI is encountered 
at steeper descent angles than the “clean” conhguration. In the -6° case, a strong 
advancing side BVI is present. In the -7.5° case, the miss distance between the rotor 
blades and vortices decreases further, and the magnitude of this BVI is increased; in 
addition, there is evidence of signihcant retreating side BVI noise directed beneath 
the helicopter. This hight condition represents the highest BVI noise levels measured 
in steady descending hight during the test program. The microphone measurement 
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Observation Angle Below Horizon, degrees 


Figure 5.18: Moving standard deviation of BVISPL at the centerline 
“virtual observer” during the 60kts, -6° flight condition. 


locations covered by the “virtual inflight array” are shown by the dotted lines on 
each hemisphere — estimated noise levels between these points are interpolated. 

Figure [5.21 and Figure 5.22 show the simulated hemispheres generated from 
the assumed analytical model using the dependent parameters identihed from the 
measured data for each case. The same general directivity trends and levels are 
observed in the modeled hemispheres as for the measured data, but the BVI directivity 
contours are much more distinct in the hemispheres generated by the analytical 
models. This is because the measured data must be time averaged in order to 
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Figure 5.19: Main rotor averaged BVISPL noise level contours for -6° 
flight path angle condition. 


compensate for variability in the flight condition during the measurement process, 
and then interpolated over the surface of the hemisphere in order to compensate 
for the relatively sparse distribution of measurements. This results in reduced 
accuracy on the underside of the hemisphere, where the measurement geometry 
changes more quickly in time and fewer periodic averages may be taken over a given 
angle. Conversely, the htted analytical models estimate noise levels in regions not 
actually measured based on physical modeling of the noise sources. As was shown 
in the wind tunnel application, more accurate estimates of the noise generated by 
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Figure 5.20: Main rotor averaged BVISPL noise level contours for -7.5° 
flight path angle condition. 


flight conditions not measured may also be made by interpolating the dependent 
modeling parameters discovered by the parameter identification process with respect 
to the non-dimensional governing parameters that classify each flight condition. For 
example, the noise levels observed at the -7.5° flight condition may be estimated using 
the modeling parameters determined from the -3° (Figure [5. 23[ ) and -9° (Figure [5. 24[ ) 
flight conditions by a linear interpolation of the dependent parameters based on the 
change in non-dimensional inflow. 

Using the standard approach of linearly interpolating SPL on an energy basis 
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Figure 5.21; Simulated BVISPL noise level contours for -6° flight path 
angle condition. 


{p%Ms) snrface of the hemisphere based on the change in flight path angle 


(e.g. RNM[TT]1 results in the hemisphere shown in Figure 5.25 It is readily apparent 
that the estimated directivity and levels produced by this “data-driven” interpolation 
approach have little relation to the BVISPL contours measured for this condition, as 
shown in Figure [5.20 The BVI hotspot is directed farther towards the retreating 
side, and the maximum BVISPL on the hemisphere is 4 dB lower than measured. 
Without a physics-based model of BVI, the simple interpolation approach is nnable 


to predict the occurrence of BVI noise along directivity angles where it did not occur 
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Figure 5.22: Simulated BVISPL noise level contours for -7.5° flight path 
angle condition. 


in the measured data. 


Figure 5.26 shows the estimated BVISPL hemisphere for the same condition 


using an interpolation of the modeling parameters inpnt to the assumed analytical 
model. By using a physics-based model, directivity patterns of BVI noise can be 
estimated which are distinct from those observed in the measnred data. The noise 
levels and directivity agree mnch more closely with the measured data. BVISPL in 
the BVI “hotspots” are about 1 dB lower than those measured. The analytical model 
also produces a more distinct directivity pattern than the hemisphere produced 
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Figure 5.23: Main rotor averaged BVISPL noise level contours for -3° 
flight path angle condition. 


directly from measurements, since the data do not need to be interpolated from sparse 
data across the surface of the hemisphere. The parameter identification scheme more 
accurately predicts noise levels and trends than simple data-driven interpolation, even 
for large changes in the governing parameters, where the conventional interpolation 
approach fails. 

Likewise, the phenomenological modeling approach introduced in this disser- 
tation can be utilized to provide more accurate extrapolation of measured data to 


other flight conditions than those measured. Figure 5.27 shows the results of using 
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Figure 5.24: Main rotor averaged BVISPL noise level contours for -9° 
flight path angle condition. 


the conventional “data-driven” approach to extrapolate the BVISPL contours for 
the -6° (Figure [5.19) and -9° (Figure 5.24) flight path angle measured and main 


rotor averaged acoustic hemispheres to the -3° (Figure 5.23) flight path angle flight 
condition. Once again, because the conventional approach lacks any information 
about the physical noise generation process, it is unable to accurately represent the 
radiated noise for conditions signihcantly different from the measured data — the 
magnitude of the estimated BVI noise is much higher than seen in the measured 
data, and the directivity much too far towards the advancing side. 
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Figure 5.25: Conventional interpolation of BVISPL noise level contours 
for -7.5° flight path angle condition. 


Figure 5.28 on the other hand, shows the results of interpolating the dependent 
modeling parameters and using the analytical noise model to estimate BVI noise levels 
at this condition. This estimate compares well with the directivity and magnitude of 


noise shown in Figure 5.23 predicting a single BVI “hotspot” in approximately the 


same direction as the measured and averaged data, with a peak level differing by 
about f dB BVISPL. In contrast to the conventional approach, the methodology 
introduced in this dissertation is able of producing a reasonable estimate of the noise 
levels at this condition by extrapolating from data collected at the -6° and -9° flight 
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Figure 5.26: Estimated BVISPL noise level contours for -7.5° flight path 
angle condition. 


path angle conditions. 
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Figure 5.27: Conventional extrapolation of BVISPL noise level contours 
for -3° flight path angle condition. 
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Figure 5.28: Estimated BVISPL noise level contours for -3° flight path 
angle condition. 
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5.3 Combined Wind Tunnel and Flight Test Model 


The ideal set of data for the modeling approach introduced in this dissertation 
would consist of high quality acoustic measurements covering a wide range of operating 
conditions — as defined by the non-dimensional governing parameters — for the specific 
full-scale rotor configuration of the modeled rotorcraft. Flight test measurements of 
helicopter noise characterize specific rotor configurations but for practical reasons 
cannot be made across the entire range of ambient and flight conditions for which 
the helicopter may be expected to operate. Wind tunnel tests allow for carefully 
controlled variation of the rotor operating condition and a high measurement quality, 
but it is impractical to conduct full-scale wind tunnel tests for each specific rotor 
configuration used on prodnction helicopters. A general-purpose rotor noise model 
should make use of both types of test data: wind tunnel data to capture the general 
parametric trends of the noise model over a wide range of operating conditions, and 
flight test data in order to tnne the model to a similar, bnt specific and fnll-scale, 
rotor configuration. 


5.3.1 Combined Modeling Approach 

A model of the OLS wind tunnel rotor was constrncted in terms of the fonr 
non-dimensional governing parameters nsing a Bayesian regnlarized 

20-nenron Artificial Nenral Network (ANN) with two hidden layers to build a 
relation between the known governing parameters and identified dependent modeling 


parameters, as described in Section 3.3.3 The dependent modeling parameters of 
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the OLS rotor ANN model were then compared against the dependent modeling 
parameters identihed for the Bell 206B3 main rotor for same non-dimensional 
operating conditions — the values of dependent parameters controlling the rotor wake 
are shown for the Bell 206B3 main rotor in Table 15.51 and for the OLS ANN derived 


model under the same conditions in Table 5.6 In general, there were no significant 


differences between the two sets of modeling parameters, with two exceptions. The 
initial vortex core size modeling parameter, geometrically scaled by rotor radius, fo, 
was signihcantly larger for the OLS rotor model than for the 206B3 model; it should 
be expected that the relative vortex core size on the OLS rotor would be larger than 
for the 206B3 because the difference in scale between the two rotors corresponds to 
a difference in Reynolds number, resulting in a difference in the non-dimensional 
vortex core size. However, the relative core size is an order of magnitude larger on 
the OLS ANN model than on the Bell 206B3 model. This difference in vortex core 


size is signihcantly larger than what would initially be expected. Figure 5.29 shows 


the variation of BVI noise at an observer ahead of the rotor as a function of vortex 
core size — from this plot it is clear that the sensitivity of BVI to vortex core size is 
much higher for values near those of the OLS ANN derived, but as the vortex core 
size decreases it’s effect on BVI noise diminishes. This is because the vortex core 
size only has a signihcant effect on the BVI process when the vortex core radius is of 
similar or greater size than the miss-distance between the blade and tip vortex. (This 


is evident from Figure 2.16 of the tangential velocity induced by a viscous vortex. ( 


Once the relevance of the vortex core parameter to BVI noise levels diminishes, the 
parameter identihcation technique of Chapter tends to select a much smaller value. 
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The other dependent parameters do not exhibit this insensitivity over the range of 



1 0 “^ 1 0 “^ 1 0 “^ 

Non-Dimensional Vortex Core Size 

Figure 5.29: Variation of peak-to-peak BVI acoustic pressure with respect 
to the non-dimensional vortex core radius. 


The second difference is a somewhat higher rate of wake contraction, modeled 
by D, for the OLS rotor model. This is a physically realistic change, as the wake is 
known to contract more quickly due to viscous effects for smaller scale rotors. [39] 


Figure 5.30 shows an estimate of the Bell 206B3 BVISPL hemisphere noise contours 


using a geometrically-scaled OLS ANN model for the -7.5° flight path angle case 
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which was previously estimated by different methods in Figures 5.20 5.22 5.25 


and 5.26 While the directivity pattern is similar to the model built from Bell 


206B3 flight test data, as shown in Figure 5.22 the magnitude of BVISPL is under 


predicted because of the relatively large vortex core size estimated from the OLS 
data. However, by adjusting the vortex core size and wake contraction parameter 
values in the OLS rotor ANN model to match those identihed for the full-scale 
206B3, a suitable full-scale model can be produced which more closely matches the 
measured full-scale data. The wake contraction term is adjusted by a bias factor, 
the best least-squares ht determined to be 0.82. As for the vortex core size, a 
linear least-squares fit would suggest all value be adjusted by dividing all values by 
a common factor of 31; although, this is an unrealistically small vortex core size. 
Instead, a more realistic factor of 5 was chosen. The output of this combined wind 
tunnel and flight test derived model is shown in Figure The combined model 


compares will with the flight test derived model, as shown in Figure 5.22, in terms 


of both levels and directivity. The result is now a model of the full-scale helicopter 
which is applicable over a wide range of operating conditions, as defined by the four 
non-dimensional governing parameters. 
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To 

2.64x10“^ 

1.98x10“^ 

1.32x10-^ 

1.30x10-^ 

c. 

8.38x10-5 

6.98x10-5 

6.34x10-5 

7.11x10-5 

7o 

1.001 

1.018 

1.003 

1.019 

7c 

-0.050 

-0.018 

-0.006 

-0.031 

7s 

0.519 

0.408 

0.001 

0.000 


Table 5.5: Dependent Parameter Values for the Bell 206B3 model. 




To 

8.00x10-3 

6.63x10-3 

5.50x10-3 

4.59x10-3 

c. 

6.34x10-3 

6.13x10-3 

6.30x10-3 

6.34x10-3 

7o 

1.020 

1.004 

1.003 

1.003 

7c 

-0.030 

-0.011 

-0.005 

-0.006 

Is 

0.522 

0.341 

0.016 

0.001 


Table 5.6: Dependent Parameter Values for the OLS ANN model. 





Figure 5.30: Geometrically-scaled OLS model BVISPL hemisphere esti- 


mate for Bell 206B3 -7.5° flight path angle condition, as in Figure 5.22 
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Figure 5.31: Combined Bell 206 - OLS model BVISPL hemisphere 
estimate for Bell 206B3 -7.5° flight path angle condition, as in Figure [5^ 
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5.3.2 Modeling the Effects of Ambient Conditions 


Existing empirical helicopter noise models — currently used in mission planning 
tools — are couched in a limited set of dimensional performance parameters, commonly 
airspeed and flight path angle. These models do not take into account the effect 
that changes in ambient conditions have on the physical process of noise generation. 
In this dissertation, phenomenological models are built in term of non-dimensional 
governing parameters. Non-dimensionalization of noise models allows them to be 
generalized to different operating conditions. 


Figure 5.32 illustrates how these non-dimensional parameters can vary with 


ambient conditions for a flight condition dehned by a constant set of dimensional 
indicated parameters — in this example a Bell 206B3 operating at a -7.5° flight path 
angle and 60 kts indicated airspeed at a variety of altitudes with ambient temperature 
and density varying as specihed by the International Standard Atmosphere (ISA) 
model. [68] Each of the four non-dimensional governing parameters influences the 


noise generation process in a different way. The effect of the individual parameter 
variations due to changes in altitude is investigated in detail in Appendix 

These variations in the non-dimensional governing parameter values can result 
in a signihcant change in the rotor acoustic state with variation in altitude. This 
effect is not accounted for in any of the empirical rotor noise modeling methods 
currently in use, all of which are developed on the basis of dimensional performance 
parameters. The effect of a change in ambient conditions can not be easily predicted 
using simple empirical corrections, since the effects of changes in the non-dimensional 
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Altitude above MSL, ft 



Figure 5.32: An example of the variation in non-dimensional governing 
parameters with altitude. 


governing parameters due to ambient conditions can be complementary or conflicting 
and depend on the specihc flight condition of the measured data, as demonstrated in 
Appendix For this particular case, the BVISPL noise level contours incorporating 
the changes in all non-dimensional parameter variations with altitude are plotted 


in Figures [5.33| and |5.34| for the 7500 ft and 15000 ft ISA cases, respectively. As 
altitude increases, BVI noise levels are shown to increase. Moreover, the directivity 
of the radiated noise is shown to shift dramatically, in the case from the retreating 


side towards the advancing side of the rotor. 
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Figure 5.33: Combined Bell 206 - OLS model BVISPL hemisphere 
estimate for Bell 206B3 -7.5° flight path angle condition at 7500 ft ISA 
altitude. 


In addition to plots of the directivity and magnitude of BVI noise, the change 
in the noise generation process with changes in the ambient conditions can be seen 
by examining the acoustic pressure time-histories at several observers for the ambient 
conditions associated with the sea level, 7500 ft and 15000 ft altitudes. The three 
selected observer locations are superimposed on the sea level BVISPL hemisphere plot 


(shown previously in Figure 5.31) in Figure 5.35 Figures 5.36, 5.37 and 5.38 show 
these pressure time-histories for the observers ahead of, toward the advancing side 
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Figure 5.34: Combined Bell 206 - OLS model BVISPL hemisphere 
estimate for Bell 206B3 -7.5° flight path angle condition at 15000 ft ISA 
altitude. 


of, and toward the retreating side of the rotor, respectively. For all three observers, 
significant changes in the pressure time-history signal are observed. The phasing of 
the BVI change with altitude in addition to the relative magnitude of the individual 


BVI pulses. For instance. Figure 5.36 shows two distinct BVI pulses at sea level, 


but as altitude increases the magnitude of the hrst BVI pulse decreases while the 


magnitude of the second BVI pulse increases. A similar trend is seen in Figure 5.37 
for the advancing side observer. For the retreating side observer, a retreating side (i.e. 
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negative) BVI pulse is seen at sea level, but weakens rapidly with altitude — at 15000 
ft, a small advancing side BVI pulse can be seen by this observer. This changes in 
acoustic pressure time-history with altitude may affect the subjective “character” of 
the BVI noise in ways not captured by the BVISPL metric. 
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Figure 5.35: Directions of observers ahead of (green), toward the advanc- 
ing side of (blue), and toward the retreating side of the rotor, superim- 


posed on the sea level hemisphere previously shown in Figure 5.31 
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Figure 5.36: Variation with altitude of the acoustic pressure time- history 
signals of one blade passage for an observer directly ahead of and 30° 
below the rotor. 
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Figure 5.37: Variation with altitude of the acoustic pressure time-history 
signals of one blade passage for an observer 30° towards the advancing 
side of and 30° below the rotor. 
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Figure 5.38: Variation with altitude of the acoustic pressure time- history 
signals of one blade passage for an observer 30° towards the retreating 
side of and 30° below the rotor. 
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The non-dimensional phenomenological models built using the methodology 
introduced in this dissertation offer a new way of estimating the changes in noise 
due to changes in the ambient condition. This is increasingly important as mission 
planning tools are being employed to estimate noise in environments considerably 
different from those where noise measurements are typically made. For example, 
while a military helicopter may be characterized at sea level conditions the noise 
models are likely to be used in hot and high conditions. Likewise, civilian helicopters 
regularly operate in all sorts of environmental extremes — for example, in the Grand 
Canyon national park alone, terrain elevations range from 2500 ft to 8000 ft above sea 
level, with temperatures ranging between 0°F and 100°F.[nn] Helicopter sightseeing 
operations are severely constrained by noise limits within national parks; the current 
practice of neglecting atmospheric effects on the helicopter noise sources is likely to 
have led to signihcantly erroneous estimates of the type and number of helicopter 
operations allowed for compliance with existing regulations. CD] 
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Chapter 6 


Conclusions 

This dissertation introduces a new non-dimensional and phenomenological 
methodology for developing analytical rotorcraft noise source models using measured 
acoustic data. In contrast, existing empirical noise source models are developed in 
terms of dimensional flight performance parameters, such as airspeed and rate of sink, 
and do not incorporate physical knowledge of the individual rotor noise generation 
mechanisms. In this new methodology, noise models are individually constructed for 
the major helicopter noise sources by applying parameter identification techniques 
to adapt assumed aeroacoustic models to measured data. The principal advantage 
of this new methodology is that it greatly enhances the applicability of measured 
rotorcraft acoustic data. This new methodology is unique in that: 

• Each of the major helicopter noise sources is considered independently 

Multiple noise sources contribute to the external noise radiation of helicopters, 
each with its own mechanisms and governing parameters. In this dissertation, 
the dominant rotor harmonic noise sources (i.e. main rotor and tail rotor 
harmonic noise) are separated from the measured data so that separate models 
can be developed for each. In addition, the rotor harmonic noise sources are 
further separated into lower harmonic loading, thickness, and blade-vortex 
interaction impulsive mechanisms. Although this dissertation focuses on main 
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rotor harmonic noise sources, especially BVI, the technique is equally applicable 
to tail rotor harmonic noise sources. 

• Noise models are developed in terms of non-dimensional governing parameters 
By developing models in a non-dimensional form, they can be generalized to 
a wide range of equivalent dimensionally dehned operating conditions. This 
concept greatly extends the applicability of the noise models; for instance, 
by applying the previonsly developed Q-SAM principle, [13] [I9] noise models 
generated from steady straight-line flight data can be utilized in order to model 
the noise radiated by gradnally accelerating or tnrning helicopters. 

• Models are built within an assumed analytical framework which incorporates 
physical knowledge of the noise generation processes 

Physically-realistic analytical models are tnned to each of the major noise 
sonrces nsing parameter identification techniqnes. Since physical knowledge 
is incorporated into the modeling framework, the generated noise models can 
be more accnrately interpolated or extrapolated to predict noise in radiation 
directions and at operating conditions for which measured data is not available. 
This phenomenological modeling approach is a significant advancement over 
the “data-driven” approaches used by existing empirical helicopter source noise 
modeling methods. 

• Models can he constructed using data from both wind tunnel and flight test 
measurements, allowing further extrapolation of the flight test data 

Both wind tunnel and flight test measurements of rotor noise have distinct ad- 
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vantages and disadvantages. Wind tunnel noise measurements can be conducted 
over a wide range of operating conditions with a high degree of consistency 
and repeatability of acoustic measurements, but full-scale measurements of 
helicopter rotors are very costly for small vehicles and impossible for large 
ones. Flight test measurements can be conducted for any rotor craft, but the 
practical range of operating conditions is limited and the flight and acous- 
tic state of the helicopter much more variable. In this dissertation, a new 
method of relating flight test data to wind tunnel data is developed, based on 
the concept of a “virtual inflight observer.” This allows the same analytical 
modeling framework to be utilized for both wind tunnel and flight test da ta. 
Moreover, this dissertation introduces a method of combining models built 
from high-quality model-scale wind tunnel measurements with models built 
from lower-quality full-scale flight test measurements of a helicopter with a 
similar rotor to produce a single model which is generalizable to a wide range 
of operating conditions, but still applicable to the full-scale vehicle. 

The new methodology developed in this dissertation was applied to measured 
wind tunnel and flight test data of two-bladed rotors. Using these data, the following 
capabilities were demonstrated: 

• Accuracy of existing methods retained 

The models constructed were accurate to within the variability of the measured 
flight test data used by existing methods. Moreover, due to the assumed ana- 
lytical structure of the modeling method, the error introduced in noise radiated 
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in directions underneath the rotor by conventional hemisphere averaging and 
interpolating schemes is limited. Lastly, by extending the applicability of high 
quality wind tunnel data to full-scale rotors, the negative effects of the inherent 
acoustic variability of helicopter flight on the modeling process can be reduced. 

• Noise estimates made in directions not measured 

Due to the assumed analytical form of the noise source models, noise can be 
estimated in directions not measured directly during testing, as was demon- 
strated using a model built from wind tunnel measurements. This capability 
is especially important when generalizing straight-line flight models to steady 
turning flight, where the helicopter will be banked, and noise radiated above the 
rotor during straight-line flight will now be radiated towards the ground. For 
practical reasons, conventional flight test measurements utilize ground-based 
microphones, making measurements of noise at or above the plane of the rotor 
impractical. 

• Noise estimates made at flight conditions not measured 

The noise models produced using the methodology introduced in this disser- 
tation are more easily interpolated or extrapolated to flight conditions not 
included in the measured data set. Wind tunnel data was used to show that 
a limited set of test data can be used to produce noise models applicable 
across a wide range of flight conditions dehned by non-dimensional governing 
parameter variations. A more limited set of flight test data was then utilized to 
demonstrate the improved generalizability of the modeling approach introduced 
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in this dissertation over conventional empirical modeling methods: the model 
constructed using the methodology developed in this dissertation was able 
to produce accurate estimates of noise at flight conditions far away from the 
data used for modeling where the conventional approach produces inaccurate 
results. 

• Noise estimates made at ambient conditions not measured 

Existing empirical models neglect the impact of ambient conditions on the 
noise generation process. Using a model constructed from wind tunnel and 
flight test data, the effect of changes in ambient conditions on the noise 
generation process is demonstrated, and found to be significant throughout a 
practical range of helicopter operating environments. Since existing empirical 
models lump all rotor noise sources together and are couched in terms of 
dimensional performance parameters, they are incapable of correcting models 
for differences between the ambient conditions of the measurement and the 
operating environment. In contrast, the non-dimensional phenomenological 
models developed in this dissertation are readily adapted to provide noise 
estimates for different ambient conditions than those measured. 

The new modeling methodology presented in this dissertation represents a 
significant advancement over the state-of-the-art in empirical rotorcraft noise source 
modeling. The applicability of measured helicopter noise data is greatly extended 
to a wide range of operating conditions while maintaining practical data collection 
requirements. The methodology introduced in this dissertation is applicable to 
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all types of rotorcraft and can be easily extended with higher hdelity modeling 
schemes for increased accnracy and generalizability as compntational resonrces 
permit. With some additional work, described in Chapter these models will 
be ready for rontine nse in mission planning tools. In addition, the techniqnes 
developed in this dissertation have application to other nses of rotorcraft external 
noise modeling. 
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Chapter 7 


Future Work 

7.1 Remaining work for Mission Planning Tool Development 

The modeling methodology developed in this dissertation is suitable for de- 
veloping helicopter noise source models for mission-planning tools; however, some 
work remains before a complete modeling tool can be developed which is suitable for 
routine use. The specific development tasks necessary to accomplish this goal are: 

• The incorporation additional physical knowledge into analytical lower harmonic 
noise source models 

In this dissertation, very simple models were developed for the lower harmonic 
noise sources of the rotor. For thickness noise, the pressure time history was 
predicted directly from the monopole displacement term in the Ffowcs Williams 
- Hawkings equation, however it is known that this formulation under-predicts 
the apparent magnitude of thickness noise seen in measurements, especially at 
higher tip Mach numbers — this is probably because the monopole thickness 
noise formulation only includes the contribution of the airfoil surface itself to the 
displacement of the medium, and neglects off-surface aerodynamic phenomena, 
such as the airfoil boundary layer and separated flow downstream of the blade 
section. The thickness noise model used in this dissertation should be expanded 
to include adjustable physics-based dependent parameters, much like the other 
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noise sources modeled in this dissertation. For example, one approach might 
be to allow the mass sinks used to model the fluid displacement caused by 
the airfoil aft of the point of maximum thickness to be redistributed past the 
trailing edge to represent separated flow. 

Likewise, lower harmonic loads were determined by assuming a simple 
spanwise triangular distribution of lower harmonic loads, and then using the 
parameter identification technique to determine the first few harmonics of the 
variation of this assumed load distribution by rotor azimuth. The advantage 
of directly specifying the lower harmonic loads in this fashion is that the 
blade motion (e.g. flapping, feathering) and full rotor trim do not need to 
be calculated in order to determine the blade loads necessary to reproduce 
the measured lower harmonic loading noise. However, because this approach 
does not incorporate much physical modeling of the lower harmonic airloads, 
the resulting lower harmonic loading noise models are not as amenable to 
accurate generalization to other operating conditions as the BVI noise model. 
A more powerful extension of the model could be developed which predicts lower 
harmonic airloads using the induced velocities determined from the wake model, 
incorporates a simple articulated rigid blade motion model, and determines the 
control input necessary to trim the rotor for each flight condition. This process 
would increase the computational costs of model evaluation, but this may be 
partially offset during the parameter identification process by a reduction in the 
number of dependent parameters needed to define the rotor state. Additional 
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refinements to the models might be made to account for the effects of active 
rotor control (e.g. higher harmonic twist/pitch control, trailing edge flaps, split 
flaps) on lower harmonic noise. These extensions would be particularly useful 
for designing helicopter operations to avoid detection, which is generally set by 
the lower frequency components of rotor harmonic noise. 


The development of a systematic approach for addressing the secondary non- 
rotor harmonic noise sources 

The approach introduced in this dissertation is valid for the dominant rotor 
harmonic noise sources, both main and tail rotor, but neglects other rotorcraft 
noise sources which can contribute to the overall external noise radiation 
under some flight conditions or rotor conhgurations. The rotor harmonic noise 


separation technique developed in Section 4^ allow the non-rotor-harmonic 
noise sources to be isolated. The contributions of these noise sources could be 
stored directly as an acoustic hemisphere, as in existing empirical methods, but 
the correct set of non-dimensional governing parameters for these noise sources 
should be identihed to ensure that the measured noise is correctly generalized. 
Alternatively, a frequency-domain parameter identification approach might be 
developed for rotor broadband noise sources and engine noise predicted for a 
given operating condition from existing empirical models. ini In either case, 
further development is required. In addition, main rotor - tail rotor interactions 
may result in signihcant external noise radiation for some rotor conhgurations 
and hight conditions; while this noise source is not harmonic with either the 
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main rotor or tail rotor, it is deterministic and might be modeled in a similar 
way as main rotor BVI noise. However, a method of isolating the contribution 
of main rotor - tail rotor interactions to the non-rotor-harmonic noise signal 
must be developed. 

• Validation of modeling method for multi-bladed rotor systems 

The measured data used in this dissertation are for two-bladed rotor systems. 
In concept, the analytical models developed in this dissertation are applicable 
to multi-bladed rotor systems. The non-dimensionalization arguments and the 
analytical acoustic model are well established for rotors with any number of 
blades. However, the addition of rotor blades leads to a more complicated wake 
structure which may not be adequately captured by the assumed Beddoes- 
derived wake model. The methodology of this dissertation should be validated 
against multi-bladed rotor data, and the assumed wake model rehned, if 
necessary, for accurate BVI noise modeling. 

• Validation of extrapolation of full-scale helicopter noise models to other ambient 
conditions than those measured 

This dissertation shows the potential of using non-dimensional phenomeno- 
logical noise models to estimate the changes in noise levels due to changes in 
ambient conditions. These results must be validated against measured flight 
test data for a helicopter operated in different ambient conditions. No suitable 
data set exists at this time, so a new flight test program must be developed 
and executed to this end. 
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Validation of the extrapolation of full-scale helicopter noise models to other 
flight conditions, such as steady turns 

Existing modeling methods are valid only for steady straight-line flight condi- 
tions. The new physics-based and non-dimensional modeling method developed 
in this dissertation is immediately applicable to slowly maneuvering flight con- 
ditions and provides a framework for future extensions to the modeling of noise 
during transient maneuvers. The Quasi-Static Acoustic Mapping (Q-SAM) 
method, [13] described in Section 


2 . 2 . 2.2 


uses momentum theory arguments to 


make an acoustic equivalence between moderate accelerations or decelerations 
and flight path angle. The Q-SAM method was later extended by the author 
to steady-turning flight conditions, allowing acoustic equivalences to be 
made for slowly maneuvering conditions terms of bank angle, acceleration, and 
flight path angle. However, under these conditions, the thrust of helicopter 
rotor can no longer be assumed to always equal the weight of the vehicle, and 
measuring the noise radiated by helicopters under these conditions is difficult. 
In addition, the Q-SAM method makes equivalencies based on the factors that 
control main rotor harmonic noise only — under some flight conditions the tail 
rotor governing parameters may vary independently of those for the tail rotor 
(e.g. tail rotor thrust and inflow may change during a turn). For these flight 
conditions where there are also signihcant levels of tail rotor noise, this main 
rotor based “equivalence” may lead to erroneous results. While the effect of 
acceleration on BVI noise has been well established, only preliminary testing of 
noise levels in steady turning flight has been conducted. |72] The extrapolation 
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of noise models to steady turning flight conditions must be more fully validated. 


• Adapt modeling method to address non-steady flight conditions, including tran- 
sient maneuvers 

For quickly maneuvering flight conditions, the quasi-static assumption of Q- 
SAM may no longer hold, and additional parameters, such as non-dimensional 
roll and pitch rates, should also be considered in any model. However, the 
method developed in this dissertation may provide a framework from which to 
develop models suitable to more rapidly maneuvering flight conditions. The 
physics-based quasi-static models built using this method may be useful in 
defining initial conditions for an aeroacoustic model with a time-varying wake 
suitable for estimating the effects of transient maneuvers. Further refinement 
could be achieved by applying parameter identification techniques to a limited 
set of data for transient flight conditions. 

7.2 Other Uses 

In addition to developing noise models for mission-planning purposes, the 
methodology developed in this dissertation could be adapted to other useful purposes. 

7.2.1 Generalization of Theoretical Noise Models 

In addition to accepting inputs from wind tunnel and flight test data, the 
modeling methodology developed in this dissertation could be utilized to construct 
lower-fidelity analytical models from the output of higher-fidelity theoretical models 
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of helicopter external noise radiation, for example combined compntational flnid 
dynamics - computational structural dynamics models. US] [ZD The development of 
accurate state-of-the-art theoretical models of helicopter noise for any given rotor 
configuration and operating condition is time consuming, but such theoretical models 
offer the only means of predicting the acoustic characteristics of new rotorcraft 
before they are constructed. Developing theoretical noise models across a wide range 
of operating conditions for a new rotor conhguration is impractical; however, the 
modeling methods developed in this dissertation could be used to generalize a limited 
set of theoretical noise models to a much wide range of operating conditions at a 
much lower computational cost. This could be useful in assessing the impact of rotor 
design changes on community and military noise metrics. 

7.2.2 Real-Time Noise Estimation 

Although the parameter identification process is time consuming, the resulting 
analytical models are lightweight and can yield estimates of external noise radiation 
in real time. Since these models are generalizable to a wide range of flight conditions 
and have no limitations in terms of directivity, they can be used to estimate noise in 
real time for arbitrary helicopter trajectories and observer locations without a priori 
knowledge of the helicopter operation. This would allow a helicopter pilot to more 
easily avoid noise sensitive areas, alert him to potential detection by adversaries, and 
provide immediate corrective feedback for maintaining a low noise flight state. 
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7.2.3 Far-Field Acoustic State Estimation 


In addition to providing a usefnl tool for pilots, real-time noise estimation 
could be extended to “close the loop” on active noise control systems. Previous 
research into active noise control for helicopters has either utilized manual sweeps of 
the control variables or closed-loop feedback control utilizing microphones located 
on the helicopter body itself. eg Open loop control is not ideal for real world noise 
control systems, since even small changes in the rotor operating condition can lead to 
large changes in noise radiation, such that any pre-programmed active noise control 
inputs are unlikely to work well under actual flight conditions. On the other hand, 
existing closed-loop approaches are wholly unsuitable as it is well known that noise 
measured on the vehicle does not correlate with external noise radiation in any direct 
fashion. The framework developed in this dissertation, however, offers a method 
of estimating the acoustic state of the helicopter in all directions. The analytical 
framework developed in this research might be extended to higher-harmonic control 
inputs — at least with respect to the lower-harmonic noise sources — allowing far-held 
noise predictions to be made. In addition, the phenomenological models developed 
in this work are applicable to both near- and far-held noise and might be used to 
relate inhight noise measurements to the overall acoustic state of the vehicle. 

7.2.4 Time-Domain Receiver Modeling 

Although most receiver models in current use are based on frequency-domain 
relations, there is recent research which indicates that these models may be inadequate 
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for representing hnman anral detection, identification, or annoyance in response 
to helicopter noise. jHj Increased interest in non-linear models of hnman hearing 
means that helicopter noise models which are capable of representing external noise 
radiation in the time-domain may soon be reqnired. Unlike any helicopter noise 
models in cnrrent nse, the modeling framework developed in this dissertation is 
capable of providing this data. 

7.2.5 Electronic Detection and Classification 

While it is still inconclusive whether the human ear can perceive phase informa- 
tion, it is readily apparent that electronic devices can. The models developed using 
the method in this dissertation could be used to train electronic detection equipment 
to recognize helicopters in a wide range of operating conditions and directions. In 
addition, using this method to extend the size of the training set will increase the 
accuracy and robustness of electronic acoustic detection and classihcation equipment. 
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Appendix A 


Theoretical Development of the Virtual Observer 

This Appendix describes the theoretical basis of the virtual inflight observer 


concept, previously developed by the author |55] and introduced in Section 4.2 of this 


dissertation. 

When a helicopter is in steady-state flight, it is creating harmonic noise that is 
associated with and necessary for powered flight. The resulting acoustic waves are 
propagated in the medium (air) and constitute the harmonic noise radiation of the 
helicopter. Motion of the helicopter has an effect on how the waves are deposited in 
the medium. Motion of the observer does not influence the sound held produced by 
the source, but does affect the acoustic pressure time-history as it is measured at 
the observer’s location. The mechanism of the Doppler effect due to motion of both 
source and observer is described in this section, using the general equations for a 
sound source in motion. I 


A.l Theory 

The general equation for the acoustic pressure observed at location x due to 
an acoustic source of strength q is given in Equation |A.1[ 


p\x,t) = 


9(9. r) 


47t|x — y \ 


(A.l) 
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Where t is the time the acoustic pressure is observed at x and r is the associated 
time at the source. 

If the source is assumed to occupy a single point in space, Xg, it may be 
expressed as follows: 

Q{y,r) = Q{T)5{y- X s{t)) (A.2) 


Substituting into Equation |A.l 


p'{x,t) = 



V 


Q(T)S(y- Xg(r)) 

4vr|a; — 


d^ydr 


(A.3) 


— oo 

By evaluating the point source at the single time of emission associated with 
the time of observation, the following expression is obtained using the retarded time 
equation: 


p{x,t) = 



Q{r)S{y - Xs{t)) 


5{t — T 




!y ATi\x-y\ ■ Oo 

— OO 

Resolving the integral in space over the compact source: 


)d'^ydT 


(A.4) 


7 Q(rM(( - r - tJM) 

— OO 

A useful property of 5-functions is given in Equation |A.6 


(A.5) 


/(r)5(^(r))dr = 


/(^) 

I 

I dr I 


(A.6) 


T* is the root of g{r) - in this case, the retarded time equation. 


g{r*)=t-r*-^-^^ = 0 
Oo 


(A.7) 
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Carrying out the differentiation of the denominator, the result is dependent 
on the Mach number of the source with respect to the medium along the radiation 


direction between the source and observer. 


dr 


-1 + 


Xi - Xsi 1 dXsi 


\x — Xs \ ao dr 


= II - M. 


(A.8) 


Substituting Equations A.6 and A.8 into Equation A. 5 results in the equation 


for the acoustic pressure generated by a moving point source. The point source is 


amplified by the Doppler amplification factor 1 1 — 


p{x,t) = 


Q{r* 


4:7t\x — Ts('T*)| |1 — Mr 


(A.9) 


Eor a source traveling along direction xi at a constant speed f/, the correct 
emission time r* can be calculated directly from a quadratic equation from the 
geometry shown in Eigure [AT| Taking the subsonic root yields Equation A. 10 

M{xi - ut) ^J{xl - my + (1 - M^){xi + xf) 


T* = t 


(A.IO) 


ao(l - M2) ao(l - M^) 

By substituting Equation |A.10| into Equation |A.9[ the acoustic pressure gener- 
ated by the moving source in the medium may be calculated for any point x at any 
time t. Motion of the observer is accounted for by dehning the location of the observer 
within the medium, x, as a function of observer time, t. It is this relation between 
source time, r*, and observer time, f, which causes the well known Doppler frequency 
shift by expanding or contracting the acoustic pressure time-history emitted by the 


source in time, as seen by the observer. (Eigure A. 2 ) 


The classical frequency-domain expression for the Doppler frequency shift may 


be obtained from the retarded time expression given by Equation A.IO For example. 
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Figure A. 2: Time dilation of pressure time-history signal of a moving 
source as seen by a stationary observer. 
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consider a simple harmonic sonrce, as described by Eqnation A. 11 Let the sonrce 
move at speed U in the same direction as the observer, moving at speed X as shown 
in Figure |A.3[ 

Q{t*) = sin(o;r*) (A. 11) 



Using Equation A. 10 the apparent frequency seen by the observer can be readily 
determined. Applying the source and observer geometry shown in Figure A. 3 to 
Equation |A.ll yields, 


M{Xt + Xo- Ut) - ^{Xt + Xo-Uty 
ao(l - M2) 

Mot + xq — Mt 1 — Mo^ I xq 

(1 - M) ^ 1-m' ao(l - M) 


(A.12) 


Substituting Equation |A.12 into A.ll yields the well known frequency- domain 
expression for the change in observed frequency due to the motions of the source 
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and observer. 


where 


Q[t) = sm(ci;- —t + uj 


Xq 


l-M ao(l-M)' 

= sm{uj't + A0) 


, 1-Mo 

^ 

l-M 


(A.13) 


However, when the direction of propagation does not lie along the direction of 
motion of both the source and the observer, the change in apparent frequency is no 
longer constant in time. Instead, the Mach numbers of the source and observer as 
measured along the direction of propagation vary over time yielding a time-varying 
change in the apparent frequency. This effect is fully accounted for in Equation |A.10 
which is valid for any motion of the observer relative to the moving source. 


A. 2 Sample Calculation 


The case of a moving point source with strength varying harmonically, as 


described in Equation A. 11 , is considered for three observer geometries. The source 


motion is similar to that of a helicopter during a level flight flyover, and is described 
in Table A.l The observer geometry is described in Table A. 2 The first observer 
geometry is that of a typical ground-based microphone, stationary with the medium. 
The second is an air-based observer which travels a fixed distance away from the 
source, but sweeps along the same elevation angles covered by the ground-based 


observer during the flyover. The third observer also travels with the source, but 
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maintains a constant elevation angle as well a distance, analogons to a stationary 
wind tunnel or inflight measurement. The ground- and air-based observer geometries 
are shown in Figure and are evaluated using Equation A. 10 The wind tunnel 
geometry is similar to the air-based geometry, except that observer travels with 
the source at a flxed angle such that the relative velocity between the observer and 
source is exactly zero. As such, the signals observed by the air-based and wind 
tunnel observers are nearly identical. 



Figure A. 4: Ground- (blue) and air- (red) based observer geometries, as 
evaluated by Equation |A.10 


Figure [A75| shows the resulting pressure time- history for the 15 Hz source at all 
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u 

170 ft/s 

Oq 

1100 ft/s 

R 

30 ft 

u 

15 & 150 Hz 

fs 

20 kHz 


Table A.l: Example calculation source motion parameters. 



Ground 

Air-Based 

Wind Tunnel 

Xi 

1000 ft 

Ut + Rcos{6g) 

Ut + i?cos(7r/4) 

X2 

0 

0 

0 

X3 

492 ft 

—Rsm{9g) 

i?sin(7r/4) 


Table A. 2; Example calculation observer geometries, where 6g is the 
elevation angle of the ground based observer, i.e. tan(6*g) = x^g/xig 
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Acoustic Pressure, Pa 


three observers for the hrst 0.2 seconds of the simulated flyover. Amplitudes of all 
three signals are approximately equal, however due to the continuously decreasing 
time delay between the moving source and the stationary ground observer, the signal 
observed on the ground is compressed in the time domain. 



Figure A. 5; Pressure time-history of 15Hz signal observed by: ground 
(blue), inflight (red), wind tunnel (black). 
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SPL, dB 



Figure A. 6: Frequency spectra of first second of 15Hz signal observed by: 
ground (blue), inflight (red), wind tunnel (black). 
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SPL, dB 



Figure A. 7: Frequency spectra of last second of 15Hz signal observed by: 
ground (blue), inflight (red), wind tunnel (black). 
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SPL, dB 



Figure A. 8: Frequency spectra of middle second of 15Hz signal observed 
by: ground (blue), inflight (red), wind tunnel (black). 
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Figures [X^|A.7[ and |A.8| show the frequency spectra of the signals as observed 
during the hrst, last, and middle second of the flyover, respectively. Dnring the 
hrst second of the flyover, the sonrce is approaching the gronnd-based observer. 
During the middle second of the flyover, the source passes overhead the ground-based 
observer. Dnring the last second the sonrce travels away from the gronnd-based 
observer. The frequency spectra are calculated using a Hann window over one second 
of data, in order to rednce spectral leakage. The air-based signal appears as a 
steady 15Hz tone in all three spectra, as does the wind tnnnel signal. However, 
the peak valne of the gronnd-based signal is shifted in freqnency, increasing as the 
source moves towards the observer and decreasing as the source moves away. More 
importantly, the frequency spectra for the ground-based signal is “smeared” across a 
broader range of freqnencies, becanse the time-delay between sonrce and observer 
varies thronghont the one second freqnency spectra window. The “smearing” is 
most severe when the sonrce passes directly overhead, since the rate of change of 
the directivity angle (i.e. slew rate) is highest. “Smearing” of the inflight spectra is 
minimal due to the slight motion between the source and observer. As expected, no 
“smearing” is observed for the stationary wind tunnel observer, producing a frequency 
spectra nearly identical to that seen by the air-based observer. 

Figures A. 9, A. 10 and A. 11 show the same frequency spectra for the 150 Hz 


source. A greater shift in frequency content is observed. 
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SPL, dB 



Figure A. 9; Frequency spectra of first second of 150Hz signal observed 
by: gronnd (blue), inflight (red), wind tunnel (black). 
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SPL, dB 



Figure A. 10: Frequency spectra of last second of 150Hz signal observed 
by: ground (blue), inflight (red), wind tunnel (black). 
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Figure A.ll: Frequency spectra of middle second of 150Hz signal observed 
by: ground (blue), inflight (red), wind tunnel (black). 
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The case of a square pulse source is also analyzed. The pressure time-history 


observed for a square pulse source with a 10% duty cycle activated at 15 Hz is shown 
in Figure A.12[ Again, the observed signal is compressed or expanded in time due to 
the change in retarded time, but the pulse shape is not distorted. 



Time, s 


Figure A. 12: Pressure time-history of square pulse signal observed by: 
ground (blue), inflight (red), wind tunnel (black). 


Figures [A.13[ [mi and | A. 15] show the hrst, last, and middle second frequency 


spectra for the square pulse source. The effect of frequency “smearing” masks the 
clear 15 Hz harmonic peaks produced by the square pulse for the ground-based 
observer, especially for the middle second where the source is flying directly over the 
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SPL, dB 


observer and the geometry between the sonrce and observer changes most rapidly. 


The tonal content of the source remains distinct for the inflight observer throughout 

the simulated flyover and is nearly identical to that observed in the wind tunnel case. 
1 50 1 1 1 1 1 1 1 1 1 

140 - 
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Frequency, Hz 


Figure A. 13: Frequency spectra of first second of square pulse signal 
observed by: ground (blue), inflight (red), wind tunnel (black). 
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150 



Frequency, Hz 


Figure A. 14: Frequency spectra of last second of square pulse signal 
observed by: ground (blue), inflight (red), wind tunnel (black). 
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SPL, dB 


150 



Figure A. 15: Frequency spectra of middle second of square pulse signal 
observed by: ground (blue), inflight (red), wind tunnel (black). 
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The “smearing” of the frequency of the source due to the Doppler effect is 
clearly evident in all of the ground-based measurements, and is especially disruptive 
for cases with impulsive harmonic noise sources — a dominant type of rotorcraft 
noise — where the various tonal frequency contributions may “smear” together. On 
the other hand, the difference between the virtual inflight and wind tunnel observers 
is negligible for typical helicopter flight test geometries. 

For this single noise source, all observers are subject to the same Doppler 
amplification, since this is determined by the Mach number of the source with respect 
to the medium. Regardless, for typical helicopter flight speeds, (M < 0.2) the effect 
of Doppler amplification due to the velocity of the helicopter through the medium is 
negligible. For example, for M = 0.2 (~130 kts at sea level) the maximum possible 
Doppler amplihcation due to overall helicopter motion is only 0.35 dB. However, 
rotor noise sources are moving with respect to the vehicle at much higher Mach 
numbers than the motion of the vehicle with respect to the medium; therefore, 
Doppler amplihcation can not be fully accounted for by considering all external 
noise radiation as originating from a single noise source affixed to the vehicle. This 
has not been accounted for in previous attempts at de-Dopplerization of rotorcraft 
noise sources |77]I7H] which have attempted to “correct” for Doppler amplihcation 
due to the motion of the vehicle as a whole, treating the rotorcraft as though it 
were a hxed-wing aircraft. Instead, the individual rotor noise source must be treated 


independently, as described in Section 2.2 
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Appendix B 


Effects of fsolated Governing Parameter Variations due to Altitude 
In this Appendix, the contribution of each of the four non-dimensional governing 
parameters to variations in helicopter external noise radiation with changes in altitude 
is independently assessed, using the combined Bell 206 - OLS model, flight condition 


and ambient condition variation from Section 5.3.2 Each of the governing parameter 


variations, shown previously in Figure 5.32 has a different effect on the magnitude 
and directivity of BVI noise. 

The thrust coefficient increases as density decreases with altitude. The thrust 
coefficient has a direct effect on the strength of the tip vortices trailed by the rotor 
blades, and hence an effect on the strength of BVI noise. Using the combined Bell 


206 - OLS ANN model developed in Section 5.3.1 , the effect of thrust on BVI noise 
can be estimated in isolation. Figure [BTT] shows the corresponding increase in BVI 
noise levels when the thrust coefficient is raised in isolation to its value at 7500 ft 


ISA altitude. Figure |B72 shows the further increase in BVI noise levels due to an 
increase in thrust coefficient to the 15000 ft value. Recall that the baseline sea level 


-7.5° case was shown for the combined Bell 206 - OLS model in Figure [VM 

In order to maintain the same indicated airspeed, dynamic pressure must be 
held constant. Since density decreases with altitude, true airspeed must be increased 
to maintain the same dynamic pressure. This corresponds to an increase in rotor 
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Figure B.l: Combined Bell 206 - OLS model BVISPL hemisphere esti- 
mate for Bell 206B3 -7.5° flight path angle condition with thrust coefficient 
at 7500 ft. 


advance ratio with altitude. The advance ratio controls the epicyclodial strnctnre 
of the rotor wake, inflnencing both the strength and directivity of BVI noise. The 
advance ratio also has a secondary effect on the lateral inflow distribntion, which 
affects the relative miss-distance of advancing versus retreating side BVI. The effect 
of changes in advance ratio on noise can be seen in isolation by estimating the noise 
levels using the combined Bell 206 - OLS ANN model. The BVISPL noise contours 


are plotted for the 7500 ft and 15000 ft advance ratios in Fignres B.3 and B.4 
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Figure B.2: Combined Bell 206 - OLS model BVISPL hemisphere esti- 
mate for Bell 206B3 -7.5° flight path angle condition with thrust coefficient 
at 15000 ft. 


respectively. The change in levels and directivity with increasing advance ratio is 
readily apparent. 

Since dynamic pressure is held constant for constant indicated airspeed, to 
hrst order the longitudinal trim of the helicopter (i.e. tip-path-plane angle of attack) 
remains unchanged with altitude. However, as altitude increases and air density 
decreases, the induced inflow through the rotor must increase in order to produce 
sufficient thrust to maintain a steady flight condition. Overall, this results in an 
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Figure B.3: Combined Bell 206 - OLS model BVISPL hemisphere esti- 
mate for Bell 206B3 -7.5° flight path angle condition with advance ratio 
at 7500 ft. 


increase in the mean inflow ratio — for flight conditions where the wake is below the 
rotor, such as the baseline -7.5° flight path angle case, this results in an increase in 
the miss-distance and a decrease in BVI noise levels. This is illustrated in Figures [B.5| 


and B.6 for the 7500 ft and f5000 ft ISA altitudes, respectively. However, if the 


wake were above the rotor for the baseline case, an increase in inflow would result in 
a decrease in miss-distance and hence an increase in BVI noise. 

Lastly, due to the decrease in ambient temperature with altitude, the speed of 
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Figure B.4: Combined Bell 206 - OLS model BVISPL hemisphere esti- 
mate for Bell 206B3 -7.5° flight path angle condition with advance ratio 
at 15000 ft. 


sound decreases, leading to an increase in the rotor tip Mach number. An increase 
in rotor hover tip Mach number results in an increase in radiated noise of all sorts, 
including BVI@. The sensitivity of noise to changes in tip Mach number increases 
as tip Mach number increases. The Bell 206B3 has a relatively low hover tip speed 
700 ft/s); for helicopters with higher tip speeds, even greater changes in noise 
levels with rotor tip Mach number can be expected, especially at high forward flight 
speeds where the advancing tip Mach number enters the transonic regime. The 
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Figure B.5: Combined Bell 206 - OLS model BVISPL hemisphere esti- 
mate for Bell 206B3 -7.5° flight path angle condition with inflow ratio at 
7500 ft. 


change in noise levels is shown in Figures B.7 and B.8 for the 7500 ft and 15000 ft 


ISA altitudes, respectively. 
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Figure B.6: Combined Bell 206 - OLS model BVISPL hemisphere esti- 
mate for Bell 206B3 -7.5° flight path angle condition with inflow ratio at 
15000 ft. 
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Figure B.7: Combined Bell 206 - OLS model BVISPL hemisphere es- 
timate for Bell 206B3 -7.5° flight path angle condition with hover tip 
Mach number at 7500 ft. 
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Figure B.8: Combined Bell 206 - OLS model BVISPL hemisphere es- 
timate for Bell 206B3 -7.5° flight path angle condition with hover tip 
Mach number at 15000 ft. 
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